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ABSTRACT 


A  hyperbaric  chamber  environment  is  simulated  using  a 
FORTRAN  computer  program.  This  program  allows  the  user  to 
input  data  files  of  varying  chamber  conditions  and  provides  a 
detailed  output  which  can  be  used  for  design  or  chamber 
operations.  Specifically,  this  program  determines  life 
support,  temperature  and  humidity  control,  and  power 
requirements.  After  developing  the  theory  and  describing  the 
program  in  detail,  complete  examples  simulating  a  hyperbaric 
chamber  at  2000  fsw,  a  submarine  lockout  chamber  at  400  fsw, 
and  a  shallow  underwater  habitat  at  50  fsw  are  demonstrated. 
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CHAPTER  I 
INTRODUCTION 

A  hyperbaric  environment  exists  in  a  facility  which  may 
be  pressurized  from  one  to  several  hundred  atmospheres. 
Edmonds  (1976)  defines  a  hyperbaric  chamber  as  an 
environmental  pressure  vessel  where  the  ambient  pressure  can 
be  increased  from  one  atmosphere  to  many  times  atmospheric 
pressure  depending  on  the  structural  capabilities  of  the 
vessel  concerned.  Hyperbaric  chambers  are  also  known  as 
recompression  chambers  and  decompression  chambers.  A 
recompression  chamber  is  used  for  compressing  a  patient  or 
diver  as  part  of  a  therapeutic  regime  or  for  training.  A 
decompression  chamber  is  used  to  decompress  a  subject  already 
exposed  to  increased  pressure.  Miller  and  Koblick  (1984) 
define  an  underwater  habitat  as  a  facility  which  acts  as  a 
survival  shelter  allowing  humans  to  survive  and  live  on  the 
ocean  floor  continuously  for  several  weeks.  Hyperbaric 
chambers,  underwater  habitats,  diving  bells,  and  submarine 
lock-out  chambers  require  unique  design  considerations  to 
account  for  the  additional  pressures  and  the  effect  on  human 
physiology  and  thermodynamic  properties. 
History 

Hyperbaric  chambers  have  had  medical  uses  for  over  3  00 
years.  Edmonds  (1976),  Miller  and  Koblick  (1984),  Shilling  et 
al.  (1976) ,  and  U.  S.  Navy  Diving  Gas  Manual  (1991)  provide  a 
chronological  history  of  hyperbaric  chambers: 


1662  -  The  first  use  of  a  hyperbaric  chamber  in  medicine  was 
made  by  a  British  physician,  Dr.  Henshaw.  This  chamber  was 
pressurized  and  depressurized  by  a  pair  of  organ  bellows. 

1830's  -  Three  Frenchmen,  Junod  (1834),  Tabarie  (1838),  and 
Pravaz  (1837)  renewed  interest  in  hyperbaric  chambers. 
Junod' s  chamber  was  a  1.6  m  diameter  sphere  capable  of 
pressures  up  to  4  atm.  He  noted  the  euphoric  effect  of 
compressed  air.  Tabarie  claimed  improvement  in  the  treatment 
of  lung  disease.  Pravaz  built  a  large  chamber  capable  of 
holding  12  patients  at  a  time. 

1850 's  -  A  rapid  rise  in  the  construction  of  hyperbaric 
chambers  was  observed.  Compressed  air  was  provided  by  steam 
engines  and  chambers  capable  of  simultaneously  treating  up  to 
50  patients  were  built. 

1860  -  The  first  hyperbaric  chamber  in  North  America  was 
built  in  Ashawa,  Canada. 

1891  -  The  first  chamber  in  the  United  States  was  built  in 
Rochester,  NY. 

1921  -  Orval  J.  Cunningham  of  Kansas  City,  MO  constructed 
the  world's  largest  hyperbaric  chamber.  It  was  5  stories 
high,  had  a  diameter  of  2  0  meters,  and  was  equipped  with 
amenities  suitable  to  a  first  class  hotel. 

Since  1920 's  -  Numerous  hyperbaric  chambers  have  been 
designed  for  treatment  of  divers.  With  the  progress  of 
hyperbaric  and  diving  medicine,  chambers  are  becoming 
increasingly  complex. 

1960 's  -  The  French  and  Americans  began  work  on  underwater 


habitats.   The  French  built  Conshelf  I  &  II  while  the  United 
States  built  Man-in-the-Sea  I  &  II  and  Sealab. 
Modern  Hyperbaric  Facilities  include: 

1)  Ocean  Simulation  Facility  (OSF)  -  located  at  the  U.  S. 
Navy  Experimental  Diving  Unit  in  Panama  City,  FL.  This 
chamber  is  one  of  the  world's  largest  man-rated  hyperbaric 
facilities  with  a  depth  capability  of  2250  fsw.  A  wide  range 
of  experiments  in  biomedical  studies,  human  factors 
verification,  and  deep  diving  systems  performance  tests  are 
conducted  in  wet  and  dry  chambers. 

2)  Naval  Medical  Research  Institute  (NMRI)  -  located  in 
Bethesda,  MD.  This  complex  is  the  U.  S.  Navy's  largest 
hyperbaric  complex  fully  devoted  to  research  and  contains 
three  dry  chambers  and  a  vertical  wetpot  rated  to  22  80  fsw  and 
two  dry  chambers  rated  to  3370  fsw.  NMRI  conducts  hyperbaric 
research  focused  on  the  biomedical  and  physiological  problems 
encountered  by  the  fleet  diver. 

3)  Duke  University's  Hyperbaric  Chamber  -  located  in 
Durham,  NC.  This  facility  is  capable  of  employing  any  desired 
gas  at  high  pressures  ranging  from  150,000  ft  of  altitude  to 
1000  fsw.  Research  is  conducted  on  acoustic  speech 
modification,  high  pressure  nervous  syndrome  and  carbon 
dioxide  effects. 

A  typical  experimental  hyperbaric  facility  is  shown  in  figure 
1.1. 

Several  factors  must  be  considered  when  designing  a 
hyperbaric  chamber.    NAVFAC  DM-39   (1982)   addresses  the 


following  design  requirements: 

1)  Pressure  vessel  design  including  material  selection, 
wall  thickness,  welding  requirements,  connections,  view  ports, 
and  doors. 

2)  Fire  protection  and  extinguishing  systems. 

3)  Electrical,  utility,  and  communication  systems. 

4)  Life  support  systems  including  mixed  gas  breathing 
systems,  carbon  dioxide  control,  and  oxygen  replenishment. 

5)  Temperature  and  humidity  control. 

6)  Power  requirements. 

The  objective  of  this  program  is  to  determine  life  support, 
temperature  and  humidity  control,  and  power  requirements. 

During  most  hyperbaric  chamber  designs,  these  factors  are 
analyzed  separately  and  compiled  into  a  finished  design.  This 
program  allows  the  user  to  input  varying  design  considerations 
and  provides  a  detailed  output  of  design  requirements.  The 
most  unique  capability  of  this  program  is  the  use  of  a 
computer  to  determine  the  thermodynamic  and  psychrometric 
properties.  Traditionally,  psychrometric  charts  have  been 
used  to  determine  dehumidif ication,  moisture  removal,  and 
ultimately  energy  requirements.  Although  specially  developed 
psychrometric  charts  are  available  for  the  higher  pressures 
required  in  a  hyperbaric  environment,  they  are  time  consuming 
and  lack  the  accuracy  of  a  computer  program.  Additionally, 
reading  values  off  these  charts  requires  non-linear 
interpolation.  Therefore,  the  necessary  relationships  were 
developed  to  replace  the  psychrometric  chart  and  simulate  a 


hyperbaric  chamber  environment. 
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1.  Divers  Swimming  and  Working  in  Wetpot  are 
Supported  by  Tenders  in  Igloo. 

2.  Igloo  is  Primary  Working,  Living  and  Resting 
Space.  "Dry"  Dives  are  also  Conducted  in 
this  Space. 

3.  Decompression  Chamber  is  used  for  Sleeping 
and  for  Treatment  of  Decompression  Sickness. 

4.  Each  Space  in  the  Facility  Can  Be  Pressurized 
Individually  or  in  Combination  With  Other 
Spaces. 


Figure  1.1  Typical  Experimental  Hyperbaric  Facility 
(Shilling  et  al.  1976) 


Objectives 

Based  upon  relevant  aspects  of  human  physiology,  this 
program  determines  the  ideal  diving-gas  mixture  along  with  an 
acceptable  range  of  diving-gas  mixtures  for  a  given  depth  or 
pressure.  Additionally,  the  required  ventilation  rate  without 
carbon  dioxide  absorbers  or  the  amount  of  absorption  material 
(sodasorb)  for  saturation  diving  is  calculated.  Various 
properties  of  the  selected  diving-gas  mixture  are  determined, 
and  by  using  thermodynamic  properties  of  the  gas,  the  heating 
and  cooling  requirements  are  calculated.   Finally,  the  energy 


requirements  for  dehumidif ication  and  reheating  to  a  desired 
temperature  are  determined  along  with  moisture  removal 
requirements . 


CHAPTER  II 
THEORY  FOR  DESIGN  OF  HYPERBARIC  ENVIRONMENTS 

Basic  Physiology 

A  basic  knowledge  of  human  physiology  is  required  to 
determine  an  acceptable  diving-gas  mixture  for  varying  depths 
or  pressures.  While  a  given  mixture  may  be  acceptable  at  one 
pressure,  it  can  cause  discomfort  or  even  be  fatal  at  other 
pressures.  The  percentage  of  oxygen  required  to  sustain  life 
decreases  as  pressure  increases  in  order  to  maintain  an 
acceptable  oxygen  partial  pressure. 

The  respiratory  system  of  the  body  is  responsible  for 
oxygen  uptake  and  carbon  dioxide  removal.  Oxygen  uptake  is 
governed  by  two  main  factors  listed  by  U.  S.  Navy  Diving  Gas 
Manual  (1971) : 

1.  Respiratory  Minute  Volume 

2 .  Effectiveness  in  which  the  lung  extracts  oxygen  from 
alveolar  gas. 

Respiratory  minute  volume,  which  is  defined  as  the  product  of 
the  respiratory  rate  (breaths/min)  and  the  tidal  volume 
(volume/breath) ,  determines  the  duration  of  exposure  of  blood 
to  alveolar  air.  The  effectiveness  in  which  the  lung  extracts 
oxygen  from  alveolar  gas  depends  upon  1)  composition  and 
volume  of  air  supplied  and  2)  partial  pressures. 

The  transfer  of  gases  takes  place  due  to  differences  in 
partial  pressures.  The  transfer  of  oxygen  occurs  with  no 
complications  if  the  oxygen  partial  pressure  is  between  0.16 


and  1.6  atm  for  short  periods  of  time.  However,  if  the  oxygen 
partial  pressure  is  below  0.16  atm,  anoxia  (oxygen  shortage) 
which  may  lead  to  fainting  and  possible  brain  damage  can 
occur.  Likewise,  if  the  oxygen  partial  pressure  exceeds  1.6 
atm,  oxygen  poisoning  which  may  lead  to  convulsions  and 
unconsciousness  may  occur.  NOAA  Diving  Manual  (1991)  states 
that  a  1.6  atm  oxygen  partial  pressure  can  only  be  tolerated 
for  up  to  4  5  minutes  and  that  the  maximum  oxygen  partial 
pressure  for  long  periods  should  not  exceed  0.6  atm.  NOAA 
Diving  Manual  (1991)  also  recommends  that  for  long  term 
exposures  in  a  hyperbaric  chamber  or  underwater  habitat  the 
oxygen  partial  pressure  of  the  diving-gas  should  be  between 
0.3  and  0.4  atm.  Figure  2.1  shows  the  percentage  of  oxygen  in 
the  diving-gas  mixture  as  a  function  of  depth  and  oxygen 
partial  pressure.  Additionally,  Figure  2.1  shows  acceptable 
mixture  ranges  for  a  given  depth  and  acceptable  depth  ranges 
for  a  given  mixture. 

Under  atmospheric  conditions,  the  partial  pressure  of 
oxygen  is  0.21  atm  since  the  percentage  of  oxygen  in  air  is 
approximately  21%.  Therefore,  at  the  surface,  divers  may 
safely  breathe  100%  oxygen  for  short  periods  (hours)  since  the 
oxygen  partial  pressure  is  only  1.0  atm,  well  below  the  1.6 
atm  upper  limit.  However,  as  depth  increases,  the  pressure 
also  increases  and  1.6  atm  can  be  guickly  exceeded  using  air 
as  the  diving-gas.  Equations  2.1  and  2.2  show  the 
relationship  between  pressure  (psia)  and  feet  of  sea  water 
(fsw) . 
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P(psia)    =  (d(fiw>+33-066)  (14.696) 

33  .  066 


d(fsw)    =  2. 25P  -  33.066  2*2 


where , 

33.066  fsw  =  1  atm 

1  atm  =  14.696  psia 
Additionally,   the  partial  pressure  of  oxygen  for  given 
conditions  can  be  determined  from  equations  2.3  and  2.4. 


PPO    =  (  d(fsw)   -n.o)  (%o,) 

2     33.066  2  23 


pp°2  -  ( P}FS^}  )  (%o2) 

14.696      2  2>4 


From  equation  2.3  it  can  be  shown  that  the  greatest  depth  in 
which  air  can  be  safely  used  as  a  diving-gas  without  exceeding 
1.6  atm  partial  pressure  is  218  fsw. 

In  addition  to  quickly  exceeding  the  allowable  amount  of 
oxygen,  the  presence  of  nitrogen  renders  air  an  undesirable 
breathing  gas  at  depth.  Under  atmospheric  conditions  the 
partial  pressure  of  nitrogen  is  low  and  no  effects  are  felt. 
However,  when  the  partial  pressure  of  nitrogen  is  sufficiently 
high,  nitrogen  narcosis  or  rapture  of  the  deep  occurs. 
Nitrogen  narcosis  presents   symptoms   similar  to   alcohol 
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intoxication  and  can  be  extremely  dangerous  while  diving  at 
great  depths.  Table  2-1  indicates  the  events  which  may  occur 
while  diving  with  high  pressure  nitrogen  at  depth  (Randall 
1992)  .  Therefore,  air  is  not  a  suitable  diving-gas  and  is  not 
generally  used  for  depths  below  150  fsw.  The  advantage  of 
using  air  when  possible  is  that  it  is  inexpensive,  and  a 
chamber  can  be  ventilated  continuously  without  the  use  of 
carbon  dioxide  absorbers  which  are  addressed  later. 


Table 

2-1 

Effects 

of  Nitrogen  Narcosis  (Randall  1992) 

100 

-  150 

fsw 

light  head,  some  euphoria 

150 

-  200 

fsw 

jovial,  some  dizziness 

200 

-  250 

fsw 

laughter  may  approach  hysteria 
simple  tasks  cannot  be  performed 
without  mistakes 

300 

fsw 

depression 

350 

fsw 

may  approach  unconsciousness 

Since  helium  is  an  inert  gas  which  eliminates  the 
narcosis  effect,  it  is  often  used  at  depths  greater  than  150 
to  2  00  fsw.  However,  helium  does  have  the  side  effects  which 
are  listed  in  the  NOAA  Diving  Manual  (1991)  .  One  of  these 
side  effects  is  speech  distortion  which  increases  with  depth. 
Speech  distortion  can  be  overcome  by  using  helium  speech 
scramblers.  A  second  problem  with  helium  is  that  its  use 
results  in  rapid  body  heat  loss  since  the  thermal  conductivity 
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of  helium  is  six  times  that  of  air.  To  counter  this  problem, 
the  diving-gas  must  be  heated  for  dives  below  800  fsw.  A 
third  problem  associated  with  helium  is  the  high  pressure 
nervous  syndrome  (HPNS)  which  causes  tremors  in  the  hands  and 
jerky  movements  of  the  limbs.  This  problem  can  be  reduced  by 
using  a  slow  rate  of  compression.  Alternatives  to  helium  are 
also  listed  in  the  NOAA  Diving  Manual  (1991)  .  These 
alternatives  include  neon  and  hydrogen.  Neon  has  the 
advantages  over  helium  of  a  lower  thermal  conductivity  and 
less  speech  distortion.  However,  neon  is  much  more  expensive 
and  creates  more  breathing  resistance.  Hydrogen  is  easier  to 
breathe  at  depth,  but  it  causes  more  speech  distortion  than 
helium,  has  a  higher  thermal  conductivity,  and  can  be 
explosive.  Other  alternatives  may  be  under  development,  but 
are  beyond  the  scope  of  this  report. 

This  program  recommends  diving-gas  mixtures  based  upon  a 
helium-oxygen  mixture.  Smaller  percentages  of  oxygen  can  be 
mixed  with  helium  to  provide  a  diving-gas  mixture  which  can 
sustain  life  without  the  risk  of  oxygen  poisoning  at  deeper 
depths . 

Gas  Mixtures 

The  ideal  diving-gas  mixture  is  determined  from  an  oxygen 
partial  pressure  of  0.35  atm.  Equation  2.5  calculates  the 
percent  oxygen  in  the  ideal  diving-gas  mixture. 
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0    35  2*5 

%0,  =    (  »■■")    x  100% 

P 

at/n 


2id«aJ 


Since  the  diving-gases  only  contain  helium  and  oxygen, 
equation  2.6  calculates  the  percent  helium  in  the  ideal 
diving-gas  mixture. 


%ifeideaJ  =  100  -  %02  .  2*6 


The  minimum  and  maximum  allowable  limits  for  oxygen  percentage 
are  calculated  by  equations  2.7  and  2.8  based  upon  oxygen 
partial  pressures  of  0.16  and  0.6  atm  respectively  as 
previously  discussed. 


2  .  7 

%o,  =  (-°-^!£)  x  ioo% 

'■atm 


2.8 

%0?   =  (-2-^.)  x  100% 

•'■atm 


If  the  maximum  allowable  oxygen  percentage  is  greater 
than  21%,  air  is  recommended  as  the  diving-gas  mixture. 
However,  in  most  cases  a  hyperbaric  chamber  operates  in  an 
environment  which  necessitates  the  use  of  a  helium-oxygen 
mixture. 

For  a  given  diving-gas  mixture,  the  pressure  (psia)  and 
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depth  (fsw)  ranges  are  calculated  using  equations  2.9  -  2.12. 


^m  "  (^)  (14.696) 


2.9 


^max  =  (-^>  (14.696) 


<4in  =  (^)  (33.066)  -  33.066 


cLax   =  <-^T>  (33.066)  -  33.066 


2.10 


2.11 


2.12 


These  ranges  are  also  shown  in  the  previous  Figure  2.1. 

Chamber  Ventilation 

Ventilation  is  required  to  remove  carbon  dioxide  and 
replenish  oxygen  to  sustain  life.  When  air  is  used,  the 
required  ventilation  rate  must  be  sufficient  to  flush  the 
chamber  to  an  acceptable  level  of  carbon  dioxide.  This  rate 
shall  be  referred  to  as  the  "free  ventilation  rate".  The  free 
ventilation  rate  is  more  than  sufficient  to  provide  enough 
oxygen  to  sustain  life.  Since  air  is  relatively  inexpensive, 
the  free  ventilation  rate  is  desirable  when  utilizing  air  and 
carbon  dioxide  scrubbers  are  not  required  to  keep  the 
environment  carbon  dioxide  free. 
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The  free  ventilation  rate  must  be  sufficient  to  maintain 
an  environment  at  or  below  acceptable  partial  pressure  levels 
of  carbon  dioxide.  As  with  oxygen,  the  partial  pressure  of 
carbon  dioxide  increases  with  depth.  For  extended  periods  of 
time,  the  maximum  allowable  partial  pressure  of  carbon  dioxide 
should  not  exceed  0.005  atm.  Figure  2.2  shows  the  maximum 
allowable  C02  partial  pressure  with  respect  to  mission 
duration.  It  is  evident  that  as  the  mission  duration 
increases,  the  smaller  the  carbon  dioxide  partial  pressure 
must  be  in  order  to  be  tolerated.  In  the  event  that  C02 
partial  pressure  exceeds  the  allowable,  hypercapnia  or  carbon 
dioxide  excess  can  occur.  Table  2-2  lists  the  four  zones  of 
hypercapnia. 


Table  2-2 

Effects  of  Hypercapnia 

Zone  I 

-  No  effect  at  low  C02  partial  pressures 

Zone  II 

-  Minor  perceptive  changes 

Zone  III 

-  Distracting  discomfort 

Zone  IV 

-  Dizziness,  stupor,  unconsciousness 

Carbon  dioxide  production  is  determined  by  equation  2.13, 


co2p  =  (jy  (o2)  (Nd)  2'13 
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where , 

C02p  =  carbon  dioxide  production  (SCFM) 

Rq  =  respiratory  quotient  (0.9) 

02c  =  oxygen  consumption  (SCFM) 

Nd  =  number  of  divers 

Oxygen  consumption  is  determined  by  the  working  conditions 

listed  in  Table  2-3. 


Table  2- 

■3 

Oxygen  Consumption  of 

Working  Conditions 

Condit 

ion 

Description 

Oo  consumption 

1 

Resting 

0.5  SLM 

/  0.0177  SCFM 

2 

Light  Work 

0.7  SLM 

/  0.0247  SCFM 

3 

Moderate  Work 

1.5  SLM 

/  0.0530  SCFM 

4 

Heavy  Work 

2.0  SLM 

/  0.0707  SCFM 

5 

Severe  Work 

3.0  SLM 

/  0.1060  SCFM 

Figure  2 . 3  shows  the  oxygen  consumption  in  Standard  Liters  per 
Minute  (SLM)  for  each  of  the  working  conditions.  Standard 
conditions  exist  at  14.696  psia  and  70°  F. 

In  order  to  determine  the  free  ventilation  rate,  a  steady 
state  condition  and  complete  mixing  are  assumed.  The 
volumetric  flow  rate  of  air  at  depth  (Vair)  in  ft3/min  is 
determined  by  equation  2.14. 
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2.14 
V        =     mRT 

aiI  "  pco2 


where, 

m  =  the  mass  flow  rate  of  C02   (lbco2/ minute) 

R  =  the  gas  constant  for  C02  (35.11  lb/lbmole-R) 

T  =  the  initial  temperature  (R) 

PC02  =  the  partial  pressure  of  C02  (psf) 

(0.005  atm  or  10.58  psf  is  recommended  for  times 
greater  than  8  hours) 
The  ventilation  rate  (Sair)  in  Standard  Cubic  Feet  per  Minute 
(SCFM)  is  calculated  by  equation  2.15. 


(V  ■   ) (P     ■  ) 
aiz  14.696 


2.15 


The  mass  flow  rate  of  C02  is  calculated  by  equation  2.16 


P  16 

(14.696psia)  ( 144 psia/psf)  (C02  ) 
m   =  E- 


(35  .11  ft-lb/ lbmole-R)   (5302?) 

where,  14.696  psia  and  530  R  are  standard  conditions. 
Equation  2.16  can  be  simplified  to  equation  2.17. 


m   =  (0.1137) {C07  ) 
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where  m  is  expressed  in  lbco2 /minute. 

Since  a  steady-state  condition  is  assumed,  it  is 
important  to  determine  how  long  it  takes  to  reach  this 
condition  at  the  beginning  of  a  mission  or  when  conditions 
change.  The  time  required  to  reach  a  95%  steady-state 
condition  can  be  derived  from  equation  2.18. 


PC02   =  <-5££)  [l-e    v<*       ] 

Vair 


2.18 


where  Vch  =  the  interior  volume  of  the  chamber  (ft3) 
Rearranging  to  solve  for  t  yields  equation  2.19. 


Vail  mRT  2 .  19 


By  substituting  in  the  values  for  PC02  (0.005  atm)  , 

R  (35.11  ft-lb/lbmole-R) ,  and  a  95%  steady-state  condition, 

equation  2.19  becomes  equation  2.20. 


f     <    V<*n-r,  (   (-95)  (.005)  (14.7)  (144)  (Vair) 

Vair  (35.11)  (mT)  n  2.20 


Carrying  out  the  multiplication  yields  equation  2.21 
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t  __   (  (-0.0167)  (V.,)  (0.2863)  (V^) 

^air  mT  2.21 


where  the  factor  0.0167  converts  t  from  minutes  to  hours. 

Carbon  Dioxide  Removal 

In  most  cases,  the  working  depth  or  pressure  requires  a 
helium-oxygen  mixture  for  reasons  previously  stated.  In  these 
cases  scrubbers  are  used  to  remove  carbon  dioxide  from  the 
chamber.  Therefore,  pure  oxygen  is  injected  but  care  must  be 
taken  to  ensure  proper  mixing  with  helium  in  the  chamber. 
This  mixing  occurs  while  the  gas  mixture  is  recirculating 
through  the  dehumidif ication  and  scrubbing  system. 

The  amount  of  pure  oxygen  required  to  sustain  the  divers 
(02r)  in  Standard  Cubic  Feet  per  Minute  (SCFM)  is  determined 
from  equation  2.22. 


a  =  (a  )  (nh)  2'22 


where  oxygen  consumption  (02c)  was  previously  determined  from 
the  diver  working  conditions.  The  amount  of  C02  produced 
(lbc02)  in  pounds  per  minute  is  calculated  by  equation  2.23. 


lbm    =  (m)   (tJ  (60)  2*23 
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where  td  is  the  mission  duration  and  the  factor  of  60  converts 
minutes  to  hours. 

NAVFAC  DM-39  (1982)  lists  three  chemicals  used  as  carbon 
dioxide  absorbents.  These  chemicals  and  their  properties  are 
shown  on  Table  2-4.  Baralyme  is  no  longer  used  in  hyperbaric 
chambers  and  Lithium  Hydroxide  is  too  expensive  and  difficult 
to  handle  for  routine  use.  Therefore,  Sodasorb  is  the 
chemical  absorbent  used  in  this  program.  The  amount  of 
Sodasorb  (LBS)  required  for  each  mission  is  calculated  from 
equation  2.24  assuming  a  50%  efficiency. 


lb  co 
s         0.245 


2.24 


Although  it  is  possible  to  obtain  higher  absorbent 
efficiencies,  canister  design,  temperature,  humidity,  and 
moisture  condensation  must  be  considered.  The  50%  efficiency 
is  a  conservative  estimate  used  for  design  purposes.  During 
actual  chamber  operations,  chemical  detection  tubes  must  be 
used  to  accurately  determine  when  to  change  carbon  dioxide 
absorbent  canisters.  The  number  of  carbon  dioxide  absorbent 
canisters  required  for  a  mission  will  depend  on  the 
capabilities  of  individual  hyperbaric  chambers. 

Diving-Gas  Mixture  Properties 

Properties  of  the  diving-gas  mixture  must  be  accurately 
determined  in  order  to  calculate  energy  and  dehumidif ication 
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Table  2-4   Characteristics  of  Carbon  Dioxide  Absorbents 
(NAVFAC  DM-39  1982) 


Absorbent 

Characteristic 

Baralyme 

Lithium 
hydroxide 

Sodasorb 

Density  (lb/ ft3) 

65.4 

28.0 

55.4 

Theoretical  C02 
absorption  (lb  C02/lb) 

0.39 

0.92 

0.49 

Useful  C02  absorption 
(based  on  50%  efficiency; 
(lb  C02/lb) 

0.195 

0.46 

0.245 

requirements  which  will  be  discussed  later.  Properties  of 
concern  include  molecular  weight,  specific  heat,  and  density 
of  the  selected  diving-gas  mixture. 

The  molecular  weight  of  a  substance  is  defined  as  the 
weight  in  pounds  of  one  pound-mole  of  the  substance.  The 
molecular  weight  of  the  diving-gas  mixtures  is  dependent  upon 
the  percentages  of  oxygen  and  helium  in  the  mixture  and 
therefore,  decreases  as  the  percent  oxygen  in  the  mixture 
decreases.  This  in  turn  decreases  the  density  which  minimizes 
the  energy  requirements.  If  the  diving-gas  is  air,  the 
molecular  weight  is  29.0  lb/ lb-mole.  Otherwise,  the  mole 
fractions  of  oxygen  and  helium  in  the  diving-gas  mixture  must 
first  be  determined.  Mole  fractions  (X)  are  determined  by 
equations  2.2  5  and  2.26.  The  sum  of  the  mole  fractions  must 
always  equal  1. 


2.25 

°2         100 
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Once  mole  fractions  are  determined,  the  molecular  weight  in 
lb/ lb-mole  can  be  calculated  from  equation  2.27. 


MW   =  [{X0){MW0)]    +  [  U„e)  (AWHe)  ]  2'27 


where , 

MW02  =  32.00  lb/ lb-mole 

MWHe  =  4.002  6  lb/ lb-mole 

The  specific  heat  at  constant  pressure  (C  ) ,  in  units  of 
BTU/ lb-mole,  of  the  diving-gas  mixture  must  be  known  for 
determining  energy  requirements.  Cp  is  defined  by  equation 
2.28. 


~.,  2,  •  2,  o 

C  =(  — ) 


where,  h  is  the  enthalpy  and  T  is  temperature. 
Given  the  mole  fractions  previously  determined,  Cp  in  BTU/ lb- 
mole  can  be  calculated  from  equation  2.29. 


cp=  r<V<V]  +  i<x*e)(cpj)  2-29 
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where , 

cp02  =  6*95  BTU/ lb-mole 

CpHe  =  4.965   BTU/ lb-mole 
When  air  is  used  as  the  diving-gas,  C   is  6.95  BTU/ lb-mole. 
Table  2-5  lists  the  molecular  weights  and  specific  heats  of 
common  helium/oxygen  mixtures  and  air. 

Table  2-5   Properties  of  Common  Diving-Gas  Mixtures 


Diving-Gas  Mixture 
(%Oo/%He) 

CP 
(BTU/ lb-mole) 

MW 
flb/lb-mole) 

1/99 

4.985 

4.283 

5/95 

5.064 

5.403 

10/90 

5.164 

6.802 

15/85 

5.263 

8.202 

20/80 

5.362 

9.602 

40/60 

5.759 

15.202 

Air 

6.950 

29.000 

Density  of  air  is  calculated  using  the  ideal  gas  equation 
of  state  since  air  will  only  be  used  down  to  150  fsw  (81.36 
psia) .  The  ideal  equation  of  state  is  shown  as  equation  2.30. 


PV =  nRuT 


2.30 
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where, 

P  =  absolute  pressure  (psf) 

V  =  volume  (cu  ft) 

n  =  number  of  moles 

Ry  =  Universal  Gas  Constant  (1544.0  ft-lb/lbmole-R) 

T  =  absolute  temperature  (R) 
Equation  2.3  0  can  be  simplified  and  rearranged  to  solve  for 
density  in  equation  2.31. 


_  (MW)  (P)  2'31 

RUT 


Density  of  helium-oxygen  mixtures  is  determined  from  the 
virial  equation  of  state.  Dymond  and  Smith  (1980)  present  the 
virial  equation  of  state  shown  as  equation  2.32. 


pym 


=  i  +  £(r)  +  C(T)     +  D(T)     + 


R  T  V  v 2  v 3  7    "\0 


where,  the  values  of  P,RU,  and  T  are  the  same  as  those  used 
above,  and 

Vm  =  molar  volume  (cu  ft/ lb-mole) 

1  =  the  first  virial  coefficient 

B(T)  =  the  second  virial  coefficient 

C(T)  =  the  third  virial  coefficient 

D(T)  =  the  fourth  virial  coefficient 
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It  is  common  practice  to  calculate  the  density  of  gases  based 
upon  only  the  first  three  virial  coefficients  since  data 
beyond  this  range  are  scarce.  Dymond  and  Smith  (1980)  state 
that  uncertainties  in  the  second  virial  coefficient  are 
relatively  small  (usually  less  than  1%)  but  are  much  larger 
for  the  third  virial  coefficient.  Additionally,  the  third 
virial  coefficient  is  much  more  temperature  dependent. 

Dymond  and  Smith  (1980)  developed  the  most  recent  table 
for  determining  the  second  and  third  virial  coefficients  for 
gas  mixtures.  Equation  2.3  3  calculates  the  second  virial 
coefficient  of  a  gas  mixture. 


2.33 

Bm^     =  BX1^    (*l)   +  Sl2(r)  <*l)  (*2>  +  B22(T)    (X2)2 


where , 

Bm(T)  =  the  second  virial  coefficient  of  the  mixture 

X-l  =  the  mole  fraction  of  the  first  pure  gas 

X2  =  the  mole  fraction  of  the  second  pure  gas 

Bll  =  the  second  virial  coefficient  of  the  first  pure  gas 

B22  =  the  second  virial  coefficient  of  the  second  pure  gas 

B12  =  the  mixed  virial  coefficient 

Equation  2.34  calculates  the  third  virial  coefficient  of  a  gas 

mixture. 
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CJT)    =  C1X1(D  (X,)3  +  3C112(D  (X1)2(X2) 

+  3C122(D  (Xx)  (X2)2   +  C222(D  (Z2)3  2.34 


where , 

Cm(T)  =  the  third  virial  coefficient  of  the  mixture 

C1X1  =  the  third  virial  coefficient  of  the  first  pure  gas 

C112  =  the  contribution  of  two  molecules  of  the  first  pure 

gas  and  one  molecule  of  the  second  pure  gas 
C122  =  tne  contribution  of  one  molecule  of  the  first  pure 
gas  and  two  molecules  of  the  second  pure  gas 
Due  to  the  absence  of  data  for  determining  the  third  virial 
coefficients  for  helium-oxygen  mixtures,   the  density  is 
calculated  based  upon  the  first  two  virial  coefficients  only. 
Tables  2-6  lists  the  virial  coefficients  for  oxygen,  helium, 
and  a  helium-oxygen  mixture  respectively  from  Dymond  and  Smith 
(1980) .   Since  this  is  the  most  recent  data  available,  it  is 
obvious  that  very  little  data  exists  to  determine  the  mixed 
virial  coefficient  for  helium-oxygen  mixtures. 

By  substituting  the  values  from  Tables  2-6  into  equation 
2.33,  the  second  virial  coefficient  in  SI  units  (cm3/g-mole) 
is  calculated  from  equation  2.35. 


B=    [(11.44)  {XHe)2]    -    [(8.8)  (XHe)  (X0)} 

-    [(16.24)  (XQ)2]  2.35 
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Table  2-6   Virial  Coefficients  of  Helium,  Oxygen,  and 
Heliox  Mixtures 
(Dymond  and  Smith  1980) 


Virial  Coeff 

icients 

TfC) 

Hel 
B 

ium 
C 

Oxygen 
B        C 

Mixture 
B        C 

-183 

10.76 

277 

-241 

-12764 

-4.4 

0 

11.87 

75.5 

-21.89 

1230 

- 

25 

11.44 

72.1 

-16.24 

1163 

- 

50 

11.37 

72.3 

-11.62 

1146 

-         .— 

B  = 

cm3/g-mo 

le       C 

=  cm6/g- 

■mole2 

In  order  to  use  this  coefficient,  it  must  be  converted  into 
the  English  units  of  ft3/ lb-mole.  This  conversion  is  made  by 
equation  2.36. 


B-    ( li!?—^)  (***£) 

28316. 85cm2  lb  2.36 

=  (0.0160)  (BSI) 


Solving  the  virial  equation  of  state  for  molar  volume  (ft3/ lb- 
mole)  yields  equation  2.37. 


2.37 

y    =  {K  +  s/K2   +AKB 

m  2 
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where  K  is  a  constant  defined  as  P^T/P. 
Finally,  density  is  calculated  by  equation  2.38 


„    MW  2-38 

P  =  IT 


where  MW  was  calculated  in  equation  2.27  and  molar  volume  was 
calculated  in  equation  2.37. 

The  U.  S.  Navy  Diving  Gas  Manual  (1971)  lists  the 
densities  of  standard  diving-gas  mixtures  at  varying 
conditions.  The  resulting  error  from  these  equations  is  less 
than  0.2%. 

Psychrometric  charts 

Psychrometric  charts  are  used  to  compute  heating  and 
cooling  requirements  along  with  dehumidif ication  and  moisture 
removal  requirements.  Standard  psychrometric  charts  can  be 
used  for  calculating  the  effects  of  heating  and  cooling  moist 
air  at  atmospheric  conditions  but  cannot  be  used  with  diving- 
gas  mixtures  at  much  greater  depths.  Although  these  charts 
are  available  for  pressures  up  to  600  psia  (approx.  1300  fsw) 
in  the  U.  S.  Navy  Diving  Gas  Manual  (1971)  ,  they  are  time 
consuming  and  lack  the  accuracy  of  this  computer  program. 
Therefore,  it  is  necessary  to  develop  relationships  which  can 
be  used  in  computer  simulation.  Additionally,  reading  values 
off  these  charts  requires  non-linear  interpolation  between  100 
psia  curves.   Figure  2.4  shows  the  psychrometric  chart  in  the 
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U.  S.  Navy  Diving  Gas  Manual  (1971)  developed  for  gas  mixtures 
at  pressures  from  100  to  600  psia. 

Given  the  dry  bulb  temperature  and  saturation  lines  at 
constant  pressure,  the  moisture  mole  fraction  (moles  water 
vapor/moles  dry  gas)  can  be  determined.  The  dry  bulb 
temperature  is  defined  by  ASHRAE  (1967)  as  the  Fahrenheit 
temperature  of  moist  gas  at  rest  with  respect  to  the 
temperature-measuring  element.  The  saturation  lines  at 
constant  pressure  represent  100%  relative  humidity  at  that 
pressure. 

The  moisture  mole  fraction  in  moles  water  vapor/mole  dry 
gas  is  determined  by  equation  2.39  if  the  pressure  of  the 
environment  and  saturation  pressure  are  known. 


X      =    P* 
ws         P   -  Ps  2.39 


where , 

X^  =  the  mole  fraction  for  water  vapor  at  saturation 

(moles  water  vapor/mole  dry  gas) 
pg  =  the  vapor  pressure  of  water  in  saturated  moist  air 

(psia) 
P  =  the  absolute  pressure  of  the  chamber  (psia) 
The  values  of  Pg  have  been  tabulated  by  ASHRAE  (1967)  for 
various  temperatures.    These  values  are  reproduced  in  Table 
2-7.   Equation  2.40  is  a  best  fit  curve  through  these  points 
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Table  2-7     Vapor  Pressure  of  Water  in  Saturated  Moist  Air  (psia) 
ASHRAE  (1967) 


Absolute  P 

ressure  (Ps) 

Temp  (F) 

Ps  (psia) 

Temp  (F) 

Ps  (psia) 

40 

0.12164 

72 

0.38856 

41 

0.12646 

73 

0.40190 

42 

0.13145 

74 

0.41564 

43 

0.13660 

75 

0.42979 

44 

0.14194 

76 

0.44435 

45 

0.14746 

77 

0.45935 

46 

0.15317 

78 

0.47478 

47 

0.15907 

79 

0.49066 

48 

0.16517 

80 

0.50701 

49 

0.17148 

81 

0.52382 

50 

0.17799 

82 

0.54112 

51 

0.18473 

83 

0.55892 

52 

0.19169 

84 

0.57722 

53 

0.19888 

85 

0.59604 

54 

0.20630 

86 

0.61540 

55 

0.21397 

87 

0.63530 

56 

0.22188 

88 

0.65575 

57 

0.23006 

89 

0.67678 

58 

0.23849 

90 

0.69838 

59 

0.24720 

91 

0.72059 

60 

0.25618 

92 

0.74340 

61 

0.26545 

93 

0.76684 

62 

0.27502 

94 

0.79091 

63 

0.28488 

95 

0.81564 

64 

0.29505 

96 

0.84103 

65 

0.30554 

97 

0.86711 

66 

0.31636 

98 

0.89388 

67 

0.32750 

99 

0.92137 

68 

0.33900 

100 

0.94959 

69 

0.35084 

70 

0.36304 

71 

0.37561 
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determined  from  PCMATLAB  (1989)  with  a  correlation  coefficient 
of  0.9999965. 


Ps  =    (1.9285  x   10"11)  T5   +    (9.86  83  x   10"10)  T4 

+  (3.3193  X   10-7)  T3  +  (1.9916  x  10"5)  T2  2.40 

+  (1.0488  x  10'3)  T  +  (2.2080  x  10"2) 


where  T  =  the  chamber  temperature  (°  F)  . 

The  U.  S.  Navy  Diving  Gas  Manual  (1971)  states  that  for 
total  mixture  pressures  up  to  100  psia,  no  correction  for  the 
effect  of  total  pressure  on  the  water-vapor  saturation 
pressure  is  used.  However,  above  100  psia,  Poynting's 
equation  is  used  as  an  estimated  correction  for  total 
pressure.  Glasstone  (1947)  derives  Poynting's  Equation  shown 
as  equation  2.41. 


Ap  m   J,(AP)(PS) 


RUT  2.41 


where  VL  is  the  specific  volume  of  the  saturated  liquid. 
Values  of  V-l  are  tabulated  by  ASHRAE  (1967)  and  are  reproduced 
in  Table  2-8  in  units  of  (ft3/lb) .  To  be  used  in  Poynting's 
Equation,  Vx  must  be  expressed  in  (ft3/mole) .  This  is  easily 
accomplished  by  multiplying  Vx  by  the  molecular  weight  of 
water  (18.0158  lb/ lb-mole)  .  Then,  the  corrected  saturation 
pressure  will  be  calculated  from  equation  2.42. 
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Table  2-8     Specific  Volume  of  Saturated  Liquid 
ASHRAE  (1967) 


vl 

Specific  Vo 

lume  (cu  ft 

/lb) 

Temp  (F) 

vi 

Temp  (F) 

Vi 

40 

0.01602 

72 

0.01606 

41 

0.01602 

73 

0.01606 

42 

0.01602 

74 

0.01606 

43 

0.01602 

75 

0.01606 

44 

0.01602 

76 

0.01606 

45 

0.01602 

77 

0.01607 

46 

0.01602 

78 

0.01607 

47 

0.01602 

79 

0.01607 

48 

0.01602 

80 

0.01607 

49 

0.01602 

81 

0.01608 

50 

0.01602 

82 

0.01608 

51 

0.01602 

83 

0.01608 

52 

0.01602 

84 

0.01608 

53 

0.01603 

85 

0.01609 

54 

0.01603 

86 

0.01609 

55 

0.01603 

87 

0.01609 

56 

0.01603 

88 

0.01610 

57 

0.01603 

89 

0.01610 

58 

0.01603 

90 

0.01610 

59 

0.01603 

91 

0.01610 

60 

0.01603 

92 

0.01611 

61 

0.01604 

93 

0.01611 

62 

0.01604 

94 

0.01611 

63 

0.01604 

95 

0.01612 

64 

0.01604 

96 

0.01612 

65 

0.01604 

97 

0.01612 

66 

0.01604 

98 

0.01612 

67 

0.01615 

99 

0.01613 

68 

0.01605 

100 

0.01613 

69 

0.01605 

70 

0.01605 

71 

0.01605 
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P    =  P     +  AP  2*42 


where  Psi  is  the  uncorrected  saturation  pressure.  The 
moisture  mole  fraction  of  water  vapor  at  saturation  (X^)  is 
then  calculated  from  equation  2.39. 

The  molar  enthalpy  of  moisture  content  at  saturation  (H) 
in  (BTU/mole  dry  gas)  is  calculated  from  equation  2.43. 


H  =    (20,000)*w  2*43 


Enthalpy  is  a  thermodynamic  property  defined  by  equation  2.44. 


2  44 

H  =  U  +  PV 


where, 

U  =  internal  energy 

P  =  absolute  pressure 

V  =  volume 
Determination  of  enthalpy  is  important  for  calculating  energy 
requirements . 

Dehumidif ication  and  Energy  Requirements 

In  order  to  determine  dehumidif ication  requirements,  the 
initial  and  desired  relative  humidities  must  be  known.  Given 
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the  initial  and  desired  relative  humidities,  the  moisture  mole 
fraction  and  molar  enthalpy  of  moisture  content  can  be 
determined  for  both  conditions  by  equations  2.45  -  2.48. 


XWi   =  (RHjiXn)  2*45 


Xw?   =  (RHd)(X„)  2'46 


Hx   =  (RHJ  (H)  2'47 


H,   =  (RHH)  (H)  2*48 


where  RJ^  and  RHd  are  the  initial  and  desired  relative 
humidities  expressed  in  fractions,  and  X^  and  H  were 
previously  determined.  Once  the  moisture  mole  fraction  at  the 
desired  relative  humidity  is  known,  the  maximum  coil 
temperature  required  for  dehumidif ication  can  be  calculated. 
Equation  2.49  calculates  the  new  saturation  pressure  (Ps) . 

U„)  (P) 


P.-    W2 


3  1  +  Xv  2.49 
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Then,  the  maximum  coil  temperature  (Tc)  in  (°F)  can  be 
determined  as  the  intersection  of  saturation  pressure  curve  at 
X^  or  from  equation  2.50. 


Tc  =    (2.4135  x  102)PSS   -    (8.4538  x  102)  Ps4 

+    (1.1635  x  103)PS3   -    (8.1910  xl02)Psz  2.50 

+    (3.5407  x  102)  Ps  +    (7 .0106) 


Equation  2.50  was  calculated  using  PCMATLAB  (1989)  and  has  a 
correlation  coefficient  of  0.999974. 

Once  the  coil  temperature  is  determined,  the  moist  gas 
mixture  enthalpies  at  the  initial  and  coil  temperatures  are 
calculated.  Equations  2.51  and  2.52  calculate  the  moist  gas 
mixture  enthalpies  in  (BTU/lb)  at  the  initial  and  coil 
temperatures  respectively. 

_   [(cp)(Td)  +/g  2*51 

a  MW 


H     m     [(Cp)(Tc)    +H2]  2'52 


c 


MW 


where  the  specific  heat  and  molecular  weight  were  previously 
determined  as  properties  of  the  diving-gas  mixture. 

The  energy  (Qd)  in  (BTU/min)  required  to  dehumidify  the 
chamber  to  the  desired  relative  humidity  is  calculated  from 
equation  2.53. 
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Qd=  PVair(Ha-Hc)  2'53 


where,  Vair  is  either  the  free  ventilation  rate  for  air 
systems  or  the  recirculation  flowrate  for  helium-oxygen 
systems.  Since  the  coil  temperature  is  generally  colder  than 
the  desired  diving-gas  mixture,  the  diving-gas  must  be 
reheated  to  a  desired  temperature.  The  change  in  enthalpy 
from  the  coil  to  a  desired  temperature  (Td)  is  calculated  by 
equation  2.54. 


[(cp)<rd-rc)]  2'54 

MW 


Then,  the  energy  required  to  reheat  the  diving-gas  (Qr)  to  the 
desired  temperature  is  calculated  by  equation  2.55. 


2.55 

Qr    =  PVailMl 


The   total   energy   (Qt)   required   by   the   system   for 
dehumidif ication  and  reheat  is  calculated  by  equation  2.56. 


2.56 

Qt   "  Od  +  Qr 


However,  energy  lost  to  heat  transfer  through  the  chamber 
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walls  depends  on  insulation  and  is  not  considered  in  this 
analysis.   It  would  be  an  additional  energy  requirement. 

Moisture  Removal  Rate 

In  all  dehumidif ication  processes,  moisture  removal  must 
be  addressed.  Due  to  the  rapid  cooling  of  the  moist  air, 
water  vapor  condenses  and  must  be  removed  from  the  chamber. 
In  order  to  determine  the  moisture  removal  rate,  the  humidity 
ratios  at  the  initial  and  desired  relative  humidities  must  be 
determined.  Equations  2.57  and  2.58  calculate  the  humidity 
ratios  in  pounds  water  vapor  per  pounds  dry  air  (lbw/lbd)  at 
the  initial  (W-J  and  desired  (W2)  relative  humidities 
respectively. 


(Xw)   (18.0158) 
K  =  -1 


'i 


MW 


0     ^R 

(Xw)   (18.0158) 
RL  =  - 


'2 


MW 


where  18.0158  is  the  molecular  weight  of  water.  The  moisture 
removal  rate  (MRR)  in  (lb/day)  can  then  calculated  by  equation 
2.59. 

MRR  =   (w1  -  w2)  Vairp(l44  0)    lb/ day  2>59 

where,  Va4r.  is  either  the  free  ventilation  rate  or  the 
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recirculation  flowrate  and  the  factor  of  1440  converts 
lb/minute  to  lb/day.  This  moisture  removal  rate  results  from 
the  dehumidif ication  process  only.  The  moisture  produced  by 
the  carbon  dioxide  absorbents  is  not  included  in  this 
analysis. 
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CHAPTER  III 
DESCRIPTION  OP  SIMULATION 

Simulation  Structure 

This  program  is  designed  to  simulate  hyperbaric  chamber 
environments  such  as  decompression  chambers,  submarine  lockout 
and  pilot  chambers,  diving  bells,  underwater  habitats,  and 
ocean  simulation  facilities  under  different  conditions.  The 
theory  and  equations  developed  in  Chapter  II  are  incorporated 
into  this  interactive  computer  program  which  is  written  in 
Fortran  77  for  a  386  personal  computer.  In  order  to  design  a 
hyperbaric  chamber  with  this  program,  the  input  listed  in 
Table  3-1  is  required.  Given  this  input,  the  program 
calculates  values  and  provides  the  output  listed  in  Table  3-2. 

A  flowchart  for  this  program  is  included  as  Figure  3.1 
and  the  source  code  is  provided  as  Appendix  B. 

Subroutine  Description 

This  program  makes  use  of  several  subroutines  to 
calculate  the  required  quantities.  A  description  of  each 
subroutine  follows  in  the  paragraphs  below. 

Subroutine  "CONVERT"  converts  depth  to  pressure.  For  all 
subsequent  calculations,  this  program  uses  absolute  pressure. 
If  the  user  chooses  to  input  a  working  depth  in  feet  of 
seawater  (fsw) ,  "CONVERT"  converts  depth  to  pressure  using 
equation  2.1.   This  value  is  returned  to  the  main  program  and 
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Table  3-1   Required  Program  Input 


Reauirement 

Units 

1) 

Working  Depth 
or  Pressure 

f  sw 

(feet 

of  seawater) 
psia 

2) 

Initial  Temperature 

o  F 

3) 

Desired  Temperature 

o  F 

4) 

Initial  Relative  Humidity 

% 

5) 

Desired  Relative  Humidity 

% 

6) 

Number  of  Divers 

N/A 

7) 

Internal  Volume  of  Chamber 

ft3 

8) 

Diver  Working  Conditions 

N/A 

9)* 

Mission  Duration 

hrs 

10) 

Recirculation  Flowrate 

ft3/min 

-  These  values  are  only  required  if  C02 
going  to  be  used  during  the  mission. 

scrubbers  are 

used  for  all  future  calculations. 

Subroutine  "GASMIX"  calculates  the  minimum  and  maximum  oxygen 
percentages  for  the  pressure  of  interest  based  upon  oxygen 
partial  pressures  of  0.16  and  0.6  atm  respectively  from 
equations  2.7  and  2.8.  If  the  maximum  acceptable  oxygen 
percentage  is  greater  than  21%,  the  ideal  diving-gas  mixture 
is  air.  All  further  calculations  are  then  based  upon  air  as 
the  diving-gas  mixture.  If  the  maximum  acceptable  oxygen 
percentage  is  less  than  21%,  the  program  determines  an 
acceptable  oxygen  percentage  using  an  IF-THEN  statement. 
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Table  3-2  Program  Output 


OutDUt 

Units 

1) 

Diving-Gas  Mixture 

%He/%02 

2) 

Oxygen  Required 

or  Ventilation  Rate 

SCFM 

3*) 

Sodasorb  Required 

lbs 

4) 

Time  to  Reach  95%  Steady-State 

hrs 

5) 

Molecular  Weight 

lb/ lb-mole 

6) 

Gas  Mixture  Density 

lb/ft3 

7) 

Maximum  Coil  Temperature 

o  F 

8) 

Dehumidif ication  Energy  Requirements 

BTU/min 

9) 

Reheat  Energy  Requirements 

BTU/min 

10) 

Moisture  Removal  Rate 

lb/day 

-  This  value  is  only  provided  if 
used  during  the  mission. 

co2 

scrubbers  are 

Subroutine  "GASRANGE"  calculates  a  pressure  and  depth  range 
based  upon  equations  2.9  -  2.12  after  an  acceptable  diving-gas 
mixture  has  been  selected.  This  information  is  useful  if  the 
same  diving-gas  mixture  is  used  for  multiple  depths.  If  the 
user  is  planning  to  operate  outside  these  limits,  a  different 
diving-gas  mixture  should  be  selected. 


Subroutine  "FREERATE"  determines  the  ventilation  rate  required 
to  keep  the  carbon  dioxide  partial  pressure  at  an  acceptable 
level  when  air  is  used  as  the  diving-gas.   First,  "FREERATE" 
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calculates  the  carbon  dioxide  production  based  upon  the 
working  conditions  and  equation  2.13.  In  order  to  maintain  a 
C02  partial  pressure  below  0.005  atm,  the  volumetric  flow 
rates  of  air  at  depth  and  at  the  surface  are  calculated  from 
equations  2.14  and  2.15. 

Subroutine  "ABSRATE"  calculates  the  amount  of  oxygen  required 
to  sustain  life,  the  amount  of  C02  produced,  and  the  amount  of 
Sodasorb  required  to  maintain  an  acceptable  environment  when 
the  diving-gas  is  a  helium-oxygen  mixture.  The  amount  of 
oxygen  required  is  calculated  by  equation  2.22  while  the 
amount  of  C02  produced  in  pounds  is  calculated  by  equation 
2.23  and  the  amount  of  Sodasorb  required  is  calculated  by 
equation  2.24.  The  user  may  wish  to  design  different  C02 
scrubbing  systems  based  upon  the  amount  of  C02  produced. 

Subroutine  "STSTATE"  calculates  the  time  required  to  reach  a 
95%  steady-state  condition.  Since  steady-state  conditions  are 
assumed,  it  is  important  to  determine  the  time  required  to 
reach  this  condition.  The  time  required  to  reach  the  95% 
steady-state  condition  is  calculated  by  equation  2.21.  This 
information  is  useful  when  beginning  a  mission  or  changing 
conditions. 

Subroutines  "MOLWT",  "SPECHEAT",  and  "DENSITY"  calculate  the 
properties  of  the  selected  diving-gas  mixture.  If  the  diving- 
gas  mixture  is  air,  the  molecular  weight  and  specific  heat  are 
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constant  values  of  29.0  lb/ lb-mole  and  6.95  BTU/ lb-mole 
respectively . 

However,  for  helium-oxygen  mixtures,  "MOLWT"  calculates 
the  molecular  weight  from  eguation  2.27  and  "SPECHEAT" 
calculates  the  specific  heat  from  equation  2.29.  This 
information  is  used  in  subsequent  calculations. 

"DENSITY"  calculates  the  density  of  the  diving-gas  for 
both  air  and  helium-oxygen  mixtures.  Since  air  is  only  used 
at  low  pressures,  "DENSITY"  determines  the  density  of  air  from 
the  ideal  equation  of  state,  equation  2.31.  When  helium- 
oxygen  mixtures  are  used,  "DENSITY"  calculates  density  from 
the  virial  equation  of  state,  equation  2.32.  The  second 
virial  coefficient  for  the  mixture  is  determined  from  equation 
2.35  and  converted  to  English  units  in  equation  2.36.  This 
value  is  then  used  to  determine  the  molar  volume  from  equation 
2.37,  and  finally  the  density  is  calculated  by  equation  2.38. 
This  value  is  used  in  subsequent  energy  calculations. 

Subroutine  "PSYCHART"  calculates  the  moisture  mole  fraction 
(X^)  and  the  enthalpy  of  moisture  content  (H)  for  the  diving- 
gas  mixture.  First,  the  vapor  pressure  of  water  in  saturated 
moist  air  (Pg)  is  calculated  from  equation  2.40.  If  the 
working  pressure  is  less  than  100  psia,  no  correction  is 
required.  However,  for  pressures  above  100  psia,  Pg  is 
corrected  using  Poynting's  equation  2.41  where  the  specific 
volume  of  the  saturated  liquid  (Vx)  is  read  off  of  Table  2-8 
using  an  IF-THEN  statement.  The  corrected  value  of  Pg  is  then 
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calculated  by  equation  2.42.  The  values  of  X^  and  H  are 
calculated  from  equations  2.39  and  2.43  respectively. 

Subroutine  "XWSH"  calculates  the  initial  and  desired  moisture 
mole  fractions  (X^  and  X^)  and  molar  enthalpies  (H-l  and  H2)  • 
These  properties  are  determined  by  the  initial  and  desired 
relative  humidities.  The  values  are  calculated  based  upon  Xwg 
and  H  by  equations  2.45  -  2.48. 

Subroutine  "COILTEMP"  calculates  the  maximum  coil  temperature 
required  to  dehumidify  the  chamber  to  the  desired  relative 
humidity.  The  moisture  mole  fraction  at  the  desired  relative 
humidity  (Xw2)  was  calculated  by  equation  2.46.  The  new 
saturation  pressure  (Pg)  is  determined  by  solving  equation 
2.39  for  Ps  which  leads  to  equation  2.49.  Given  Pg/  the  coil 
temperature  (Tc)  is  calculated  by  equation  2.50. 

Subroutine  "MGMENTH"  calculates  the  moist  gas  mixture 
enthalpies  entering  the  cooling  coil  (Ha)  and  leaving  the 
cooling  coil  (Hc)  by  equations  2.51  and  2.52  respectively.  In 
this  program,  the  temperature  entering  the  cooling  coil  is 
assumed  to  be  the  initial  temperature  (TL)  and  the  temperature 
leaving  the  cooling  coil  is  assumed  to  be  the  maximum  coil 
temperature  (Tc) .  These  values  may  differ  in  an  actual 
cooling  system. 

Subroutine   "ENERGY"   calculates  the  energy  required  for 
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dehumidif ication  and  reheat  to  the  desired  temperature.  The 
energy  required  for  dehumidif ication  (Qd)  is  calculated  by 
equation  2.53.  If  the  diving-gas  is  air,  the  free  ventilation 
rate  (Vair)  calculated  by  "FREERATE"  is  used.  For  helium- 
oxygen  mixtures,  the  recirculation  flowrate  input  by  the  user 
is  used  as  Vair.  In  order  to  determine  the  energy  required  to 
reheat  to  a  desired  temperature  (Qr) ,  the  change  in  enthalpy 
from  the  coil  temperature  to  a  desired  temperature  must  be 
calculated  by  equation  2.54.  Then,  Qr  is  calculated  in  a 
similar  manner  as  Qd  by  equation  2.55.  The  total  energy 
required  (Qt)  is  calculated  by  equation  2.56. 

Subroutine  "MOISTREM"  calculates  the  moisture  removal  rate 
required  for  the  above  dehumidif ication.  First,  the  humidity 
ratios  at  the  initial  (W-^)  and  desired  (W2)  relative 
humidities  are  calculated  by  equations  2.57  and  2.58.  Then, 
the  moisture  removal  rate  is  calculated  by  equation  2.59  where 
again,  Vair  is  the  free  ventilation  rate  for  air  or  the 
recirculation  flowrate  for  helium-oxygen  mixtures.  The 
calculations  performed  by  these  subroutines  are  incorporated 
into  the  main  program,  and  the  results  are  displayed  to  the 
user  as  program  "output".  Three  examples  are  discussed  in 
detail  in  the  next  chapter. 
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CHAPTER  IV 
APPLICATION  OF  SIMULATION 

This  program  can  be  used  not  only  to  simulate  a 
hyperbaric  chamber  but  also  a  submarine  diver  lock-out  chamber 
or  an  underwater  habitat.  The  following  examples  demonstrate 
this  program's  usefulness  for  determining  the  requirements  of 
these  cases. 

Example  1:   Hyperbaric  Chamber 

This  example  simulates  a  hyperbaric  chamber  used  for  a 
deep  dive.  As  a  model,  the  U.  S.  Navy's  Ocean  Simulation 
Facility  at  the  Naval  Coastal  Systems  Center  in  Panama  City, 
Fl  was  be  used.  This  facility  as  shown  in  Figure  4.1  is 
capable  of  pressurization  down  to  2250  fsw  and  is  described  in 
detail  by  Shilling  et  al.  (1976)  and  U.  S.  Navy  Diving  Manual 
(1991) .  The  input  conditions  for  this  example  are  given  in 
Table  4-1.   The  calculation  steps  are  demonstrated  below: 

1)   Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.   By  equation  2.1, 


(2000  t   33.066)  }  (14>696) 
33.066 

P  =   903  .6  psia 
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T 

8  FT 
DIA 

1 


30  FT 


Figure  4.1  Ocean  Simulation  Facility 
(Shilling  et  al.  1976) 


2)   The  ideal  gas  mixture  for  this  depth  was  calculated 
by  equations  2.5  and  2.6. 


%0>  °  ( ;  "pa"  j ' (100%) 

14.696 
%0o  =  0.57% 


%He  =  100  -  0.57 
%He  =  99.43% 
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Table  4-1   Hyperbaric  Chamber  Simulation  Input 


Depth  (fsw) : 

2000 

Initial/Desired  Temperatures  (°  F)  : 

90/80 

Initial/Desired  Relative  Humidities  (%) : 

90/50 

Number  of  Divers: 

6 

Internal  Chamber  Volume  (ft3) : 

3300 

Recirculation  Flowrate  (cfm) : 

100 

Mission  Duration  (hrs) : 

240 

Diver  Working  Conditions:           Light  Work 

Therefore,  the  ideal  gas  mixture  was  determined  to  be 
0.57%  02  and  99.43%  He. 

3)    The  minimum  and  maximum  oxygen  percentages  were 
calculated  by  equations  2.7  and  2.8. 


%o2   =  ( B°'16t   )  (100%) 

2»i"     ,  903  .6  v 

14.696 


%09   =  0.26% 


%o2  =  ( — °;*   )  (100%) 

2~*     ,  903.6  v 
14.696 


%09   =0.9  8% 
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Since  only  mixtures  within  these  levels  are  acceptable, 
air  could  not  be  used  and  a  0.5%  02  and  99.5%  He  diving- 
gas  mixture  was  selected.  Therefore,  carbon  dioxide 
scrubbers  were  used  and  a  computation  of  free  ventilation 
was  not  needed. 

4)   The  pressure  and  depth  ranges  for  the  gas  mixture 
were  calculated  from  eguations  2.9  -  2.12. 


*»m  "  (-^)  (14.696) 


*»in  =  47  0.3  psia 


^max  "  (4r^>  (14.696) 

max       q  _  5 

^max  ■  1763.5  psia 


^in  =  t(^)  (33.066)  -  33.066] 


*mim   =    1025-°  fSW 


^nax  "  [(^)  (33.066)  -  33.066] 
d^  =   3934.9  fSW 


5)  Since  a  helium-oxygen  mixture  was  used  with  carbon 
dioxide  scrubbers,  the  amount  of  pure  oxygen  injected 
into  the  chamber  was  calculated  by  equation  2.22. 
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02c   =0.7    SLM  or   0.0247    SCFM   for    light  work 


02     =    (O.0247) (6) 


02z  =  0.14  82  SCFM 


6)  The  amount  of  carbon  dioxide  produced  and  the 
scrubber  requirements  were  calculated  by  equations  2.17 
and  2.2  3  -  2.24. 

m   =  (0.1137)  (0.9)  (0.1482) 
m  =  0.01517  lbco  /min 


lbC02  =    (0.0157)  (240)  (60) 
lbC02   =  218.4  lb 


7,     ,  218. 4N 
lbs~    ("0T245) 

lbs  =   891.4  lb 


7)  The  time  required  to  reach  95%  steady-state  was 
determined  from  equation  2.21  where  Vair  was  determined 
from  2.14. 
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v        =     [(0.01517)  (35.11)  (549.67)] 
air  10.58 

Vair  =  27.7  ft3/min 


t  =    (-0.0167)  (H£0)ln[l-      (0-2863)  (27.7)      j 
27 .7  (0.0152)  (549  .67) 

t  =  5.96  his 


8)   Next,  the  diving-gas  properties  were  determined. 
Molecular  weight  was  calculated  from  equations  2.2  5-2.2  7 


Xn    =   -°-^  =  0.005 
°2  100 


XMtk  =  21^1   =  0.995 
He  100 


MW  =    [(0.005)  (32.00)  +  (0 . 995)  (4  .  0026) ] 
MW  =  4.143  lb  I  lb-mole 


Specific  heat  was  calculated  from  equation  2.29. 
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Cp  =    [(0.005)  (6.95)  +  (0.995)  (4.965)] 
Cp  =  4.97  5  BTU/ lb-mole 


Density  was  calculated  from  equations  2.35  -  2.38 


B  =    [(11.44)  (0.995)2   -    (8.8)  (0.995)  (0.005) 
-    (16.24) (0.005)2] 

B  =  11.2187  cm3/g-mole 


B  =    (0.0160) (11.2817) 
B  =  0.1805  ft3/ lb-mole 


K  _     (1544)  (549.67) 
(903.6)  (144) 

K  =  6.5226 


„  6.5226    +  y/6.52262   +  4(6.5226)  (0.1805 

m  2 

Vm  =  0.6185  ft3 /lb-mole 


4.143 
P     '    6.6984 

p   =  0.6185   lb/ft3 
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9)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.2. 

Ps  =  0.69838  psia   (from  Table  2-7) 

VL  =  0.01610  cu  ft/lb   (from  Table  2-8) 

V-l  =  0.2880  cu  ft/mole 


AP  [(0.69838)  (0.2880)  (903.6)  (144)] 


[  (1544)  (549.67) ] 
APS  =  0.0308  psia 


Ps  =  0.69838   +  0.0308 


P„  =  0.7292  psia 


_  0.7297 

ws  "    903.6    -  0.7297 

Xws  =  0.000  81  moles  water  vapor/ mole  dry  gas 


(X^s   is  shown  as  point  A  on  Figure  4.2) 


Xwl   =    (0.9)  (0.00081) 

Xwl  =  0.00072  moles  water  vapor  /mole  dry  gas 
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Moisture  Mole  Fraction 
(moles  water  vapor/mole  dry  gas) 
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(X^  is  shown  as  point  B  on  Figure  4.2) 

Xw2  =    (0.5)  (0.00081) 

Xw2  =  0.00041  moles  water  vapor/ mole  dry  gas 


(X^2    is  shown  as  point  D  on  Figure  4.2) 


H  =    (20,000)  (0.00081) 

H  =  16  .2  BTU/mole  dry  gas 


Hx  =    (0.9) (16.2) 

Hx  =  14.6  BTU/mole  dry  gas 


H2  =    (0.5) (16.2) 

H2  =   8.1  BTU/mole  dry  gas 


10)  The  maximum  coil  temperature  was  determined  from  the 
psychrometric  chart  as  the  intersection  of  X^  and  the 
100%  saturation  curve.   This  point  is  at  approximately 
70°  F  and  is  shown  as  point  C  on  Figure  4.2. 
Tc  =  70°  F 
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11)    The  energy  requirements  were  calculated  from 
equations  2.51  -  2.56. 


H    =     [4.975) (90)    +  14.6] 
a   '  4.143 


Ha  =  111.6  BTU/lb 


H    =     [(4.975)  (70)    +  8.1] 
c  4.143 

Hc  =  86  .  0  BTU/lb 


Qd  =    (0.6185)  (100)  (111.6    -  86.0) 
Qd  =  1583.4  BTU/min 


Ah  =     [4.975)  (80   -70)] 
4.143 

Ail  =  12.0  BTU/lb 


Qr  =    (0.6185)  (100)  (12.0 
Qr  =  742.2  BTU/min 


Qt  =  1583.4    +  742.2 
Qt  =  2325.6  BTU/min 
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12)    The  moisture  removal  rate  was  calculated  from 
equations  2.57  -  2.59. 


K  =  (0.00072)  (  18-0159  ) 
1  4.143 


W1  =   0.00313  lb  J  lb  d 


W2  =   (0.00041)  (1|4tP) 
z  4.143 


W2  =   0.0017  8  lb  J  lb 


d 


MRR  =    (0.00313  -  0.00178) (100) (0.6185) (1440) 
MRR   =  120.1  lb  I day 


The  actual  program  printout  is  provided  as  Table  4-2.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  the 
example. 
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Table  4-2   Output  for  Hyperbaric  Chamber  Simulation 


HYPERBARIC  CHAMBER  SIMULATION 
Input  Conditions: 
Pressure  (psia) : 

Initial/Desired  Temperature  (°  F)  : 
Initial/Desired  Relative  Humidity  (%) : 
#  of  Divers/Working  Conditions: 
Chamber  Volume  (ft3) : 
Flow  Rate  (ft3/min) : 
Mission  Duration  (hrs) : 

Output  Conditions: 

Diving-Gas  Mixture  (%He/%02)  : 

Oxygen  Required  (SCFM) : 

C02  Produced  (lbs) : 

Sodasorb  Required  (lbs) : 

Time  to  95%  Steady-State  (hrs) : 

Molecular  Weight  (lb/ lb-mole) : 

Gas  Mixture  Density  (lb/ ft3)  : 

Maximum  Coil  Temperature  (°  F)  : 

Dehumidif ication  Energy  (BTU/min) : 

Reheat  Energy  (BTU/min) : 

Total  Energy  Requirement  (BTU/min) : 

Moisture  Removal  Rate  (lb/day) : 


903.6 
90.0/80.0 
90.0/50.0 
6.0/2 
3300.0 
100.0 
240.0 


99.5/0.5 
0.148 
218.7 
892.7 
6.0 
4.143 
0.6184 

70.2 

1563.5 

724.3 

2287.8 

125.2 
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Example  2;   Submarine  Diver  Lock-out  Chamber 

This  example  simulates  a  diver  lock-out  chamber  aboard  a 
submarine  used  to  decompress  divers  after  a  mission.   The 
input  conditions  for  this  example  are  given  in  Table  4-3. 
Table  4-3   Diver  Lock-out  Chamber  Simulation  Input 


Depth  (fsw) : 

400 

Initial/Desired  Temperatures  (° 

F) 

85/78 

Initial/Desired  Relative  Humidities 

(%) :   90/60 

Number  of  Divers: 

2 

Internal  Chamber  Volume  (ft3) : 

100 

Recirculation  Flowrate  (cfm) : 

10 

Mission  Duration  (hrs) : 

4 

Diver  Working  Conditions: 

Resting 

The  calculation  steps  are  demonstrated  below: 

1)   Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.   By  equation  2.1, 


P-  (  (400  +  33.066)  )(14.696) 

33.066 

P  =  192  .5  psia 


2)      The  ideal  gas  mixture  for  this  depth  was  calculated 
by  equations  2.5  and  2.6. 
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%02  =  ( 2^3J )  (100%) 

(  192.5  v 

14.696 


%02  =  2.67% 


%He   =  100  -  2.67 
%He   =  97  .33% 


Therefore,  the  ideal  gas  mixture  was  determined  to  be 
2.67%  02  and  97.33%  He. 

3)    The  minimum  and  maximum  oxygen  percentages  were 
calculated  by  equations  2.7  and  2.8. 


%o2.  =  ( — °-16   )  (ioo%; 

2min     ,  192  .5  v 
14.696 


%o,   =1.2% 


%°2  =  ( — ,°;6C  (ioo%) 

2**x  ,    192  .5  v 

14.696 


%07   =  4.58% 

■'max 


Since  only  mixtures  within  these  levels  air  acceptable, 
air  could  not  be  used  and  a  3%  02  and  97%  He  mixture  was 
selected.  Therefore,  carbon  dioxide  scrubbers  were  used 
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and  a  computation  of  free  ventilation  rate  was  not 

needed. 

4)       The    pressure    and   depth    ranges    for    the    gas    mixture 

were  calculated   from  equations   2.9   -   2.12. 


Pmin=    (^y)  (14-696) 


^mm  =  78.4  psia 


^max=     <-~y>   (14-696) 

^max  =  29  3.9  psia 


<3nin  =    [(yy)  (33.066)    -  33.066] 
4nin    =    143.3    fsW 


^nax  =    [(-y^)  (33.066)    -  33.066] 

<3nax    =    628.2    fSW 


5)   Since  a  helium-oxygen  mixture  was  used  with  carbon 
dioxide  scrubbers,  the  amount  of  pure  oxygen  injected 
into  the  chamber  was  calculated  by  equation  2.22. 
02c  =0.5  SLM  or  0.0177  SCFM 
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02     =  (0.0177)  (2) 


02     =0.03  54  SCFM 


6)  The  amount  of  carbon  dioxide  produced  and  the 
scrubber  requirements  were  calculated  by  equations  2.17 
and  2.2  3  -  2.24. 

m  =    (0.1137) (0.9) (0.0354) 
m   =  0.0036  lbco  /min 


lbC02   =  (0.0036)  (4)  (60) 
lbco    =   0.86  lb 


lbs=    (  °-86  ) 
s  0.245 

lbs  =  3 . 5  lb 


7)  The  time  required  to  reach  95%  steady-state  was 
determined  from  equation  2.21  where  Vair  was  determined 
from  2.14. 


v  [  (0.0036)  (35.11)  (544.67) ] 

air   '  10.58 


Vair  =  6  .51  ft3/min 


66 


t=  (-Q.Q167)  (^QQ-)ln[l-   (0-2863)  (6.51)  . 
6.51        (0.0036) (544.67) 

t   =  0.77  his 


8)   Next,  the  diving-gas  properties  were  determined. 
Molecular  weight  was  calculated  from  equations  2.2  5 
2.27. 


3  0 
Xn    =  ^-^  =0.03 
°2    100 


97  0 
X„e   =  Z±L}L   =0.97 

He  100 


MW  =    [  (0.03)  (32.00)  +  (0.97)  (4.0026) ] 
MW  =   4.843  lb  I  lb-mole 


Specific  heat  was  calculated  from  equation  2.29 


Cp  =    [  (0.03)  (6  .95)  +  (0.97)  (4.965) ] 
Cp  =   5.025  BTU/ lb-mole 


Density  was  calculated  from  equations  2.35  -  2.38 
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B  =    [(11.44)  (0.97)2   -    (8.8)  (0.97)  (0.03) 
-    (16.24) (0.03)2] 

B  =  10.4932  cm3  / g-mole 


B  =    (0.0160) (10.492) 
B  =  0.1679  ft3/ lb-mole 


K  _     (1544) (544.67) 
(192.5)  (144) 

K  =  30.34 


30.34  +  y/30.542  +  4(30.54)  (0.1679) 
m  2 

Vm   =  30.51  ft3/ lb-mole 


4.843 

P  =  

K    30.51 

p  =  0.1587  lb/ ft3 

9)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.3. 

Pg  =  0.59604  psia  (from  Table  2-7) 

V1   =   0.01609  cu  ft/lb  (from  Table  2-8) 

V1   =  0.2896  cu  ft/mole 
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Moisture  Mole  Fraction 
(moles  water  vapor/mole  dry  gas) 
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Ap     _     [(0.59604)  (0.2896)  (192.5)  (144)] 
s  '  [  (1544)  (544.67) ] 

APS  =  0.00569  psia 


Ps  =  0.59604    +  0.00569 


Ps  =  0 .6  017  3  psia 


X      =  0.60173 

ws  '     192.5   -  0.60173 


X      =  0.00314  moles  water  vapor  /mole  dry  gas 


(X^  is  shown  as  point  A  on  Figure  4.3) 

Xwl   =    (0.9)  (0.00314) 

Xwl  =   0.0  0282  moles  water  vapor /mole  dry  gas 


(X^  is  shown  as  point  B  on  Figure  4.3) 

Xw2  =    (0.6)  (0.00314) 

Xw2  =  0.00188  moles  water  vapor  /mole  dry  gas 
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(X^  is  shown  as  point  D  on  Figure  4.3) 


H  =    (20,000) (0.00314) 

H  =  6  2.71  BTU/mole  dry  gas 


Hx  =    (0.9)  (62.71) 

Hx  =  56.44  BTU/mole  dry  gas 


H2   =    (0.6)  (62.71) 

H2   =  37  .63  BTU/mole  dry  gas 


10)  The  maximum  coil  temperature  was  determined  from  the 
psychrometric  chart  as  the  intersection  of  X^  and  the 
100%  saturation  curve.  This  point  is  at  approximately 
70°  F  and  is  shown  as  point  C  on  Figure  4.3. 

Tc  =  70°  F 

11)  The  energy  requirements  were  calculated  from 
equations  2.51  -  2.56. 


„  =  [(5.025)  (85)  +  56.44] 
a  "  4.843 


H.  =   99 .85  BTU/lb 
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„  [  (5.025)  (70)    -*■  37  .63] 

4.843 

Hr  =  80.40  BTU/lb 


Qd  =    (0.1587)  (10)  (99.85   -  80.40) 
Qd  =  3  0.9  BTU/min 


Ah  =     [(5.025)  (78   -  70)] 
4.843 

Ah  =  8.3  BTU/lb 


Qr  =    (0.1587)  (10)  (8.3) 
Qr  =  13.2  BTU/min 


Qt  =  30.9   +  13 .2 
Qt  =  44.1  BTU/min 


12)         The    moisture     removal     rate     was     calculated     from 
equations    2.57    -    2.59. 


W.    =    (0.00282)  (  18.0158) 
1  4.843 


W1  =  0.01049   lbjlbd 
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W2   =    (0.00188)  (  1!'^o8) 

z  4  .  843 

W2   =  0.00699  lb  J lbd 


MRR   =  (0.01049  -  0.00699) (10) (0.1587) (1440) 
MRR  =  8.0  lb/ day 


The  actual  program  printout  is  provided  as  Table  4-4.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  the 
example. 
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Table  4-4   Output  for  Diver  Lock-out  Chamber  Simulation 


DIVER  LOCK-OUT  CHAMBER  SIMULATION 

Input  Conditions: 

Pressure  (psia) : 

Initial/Desired  Temperature  (°  F)  : 

Initial/Desired  Relative  Humidity  (%) : 

#  of  Divers/Working  Conditions: 

Chamber  Volume  (ft3) : 

Flow  Rate  (ft3/min) : 

Mission  Duration  (hrs) : 

Output  Conditions: 

Diving-Gas  Moisture  (%He/%02)  : 

Oxygen  Required  (SCFM) : 

C02  Produced  (lbs) : 

Sodasorb  Required  (lbs) : 

Time  to  95%  Steady-State  (hrs) : 

Molecular  Weight  (lb/ lb-mole) : 

Gas  Mixture  Density  (lb/ ft3) : 

Maximum  Coil  Temperature  (°  F)  : 

Dehumidif ication  Energy  (BTU/min) : 

Reheat  Energy  (BTU/min) : 

Total  Energy  Requirement  (BTU/min) : 

Moisture  Removal  Rate  (lb/day) : 
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192.5 

.0/78.0 

.0/60.0 

2.0/1 

100.0 

10.0 

4.0 

.0/3.0 
0.035 
0.9 
3.5 
0.8 
4.843 
0.1587 
70.0 
30.9 
13.2 
44.1 
8.0 
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Example  3;   Underwater  Habitat 

This  example  simulates  a  shallow  underwater  habitat.  The 
French  Conshelf  II  is  used  as  a  model.  Conshelf  II,  shown  as 
Figure  4.4,  was  a  shallow  water  habitat  operating  in  the  Red 
Sea  in  September  of  1963.  The  breathing  gas  was  air  and  the 
mission  duration  was  1  month.  Miller  and  Koblick  (1984) 
describe  this  mission  in  detail.  The  input  conditions  for 
this  example  are  given  in  Table  4-5.  The  calculation  steps 
are  demonstrated  below: 

Table  4-5   Underwater  Habitat  Simulation  Input 


Depth  (fsw) :  50 

Initial/Desired  Temperatures  (°  F) :  95/80 
Initial/Desired  Relative  Humidities  (%) :  95/50 
Number  of  Divers:  4 

Internal  Chamber  Volume  (ft3) :  2000 

Diver  Working  Conditions:  Moderate  Work 


1)   Since  a  working  depth  was  input,  it  was  converted  to 
a  pressure.   By  equation  2.1, 


P=  (  (5Q  +33-066))  (14.696) 
33.066 


P  =  36 .9  psia 
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7     9 


Figure  4.4   Conshelf  II  Habitat,  Starfish  House 
(Miller  and  Koblick  1984) 


2)  The  ideal  gas  mixture  for  this  depth  is  calculated  by 
equations  2.5  and  2.6. 


%o2  =  ( °;3*   (100%) 

(  36-9  ) 

14.696 
%0,  =  13.9% 


%He  =   100  -  13.9 
%He  =   86  . 1% 
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3)    The  minimum  and  maximum  oxygen  percentages  were 
calculated  by  equations  2.7  and  2.8. 


%0,   =  ( 2jJ^ (100%) 

2»i„       36.9  v 

14.696 


%07       =6.37% 

zmin 


%02   =  ( 1^ (100%) 

2*»x  ,     36.9  v 

14.696 


%0.   =  23.9% 


Since  any  mixture  within  these  levels  is  acceptable,  air 
(21%  02)  was  selected  as  the  diving-gas.  Therefore,  a 
free  ventilation  rate  was  used  to  keep  the  chamber  free 
of  carbon  dioxide  and  C02  absorbers  were  not  required. 

4)   The  free  ventilation  rate  at  depth  and  from  the 
surface  was  calculated  from  equations  2.14  -  2.17. 
02c  =  1.5  SLM  or  0.053  0  SCFM  for  moderate  work 

CO-     =  (0.9) (0.0530) (4) 
CO?    =0.1908  liVmin 
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m   =  (0.1137)  (0.1908) 
m  =  0.0217  lbco  /min 


„   =  [(0.0217)  (35.11)  (529.67) ] 
air  "  10.58 


Vair   =  39.9  ftVmin 


s        _  [(39.9)  (36.9)  (529.67)] 
air  "     [(14.696)  (554.67)  ] 

Sair   =  95.8  SCFM 


5)    The  time  required  to  reach  a  95%  steady-state 
conditions  was  determined  from  equation  2.21. 


t  =  (-0.0167)  (2000)ln[1    (0.2863)  (39.9)  j 
39.9        (0.0217)  (554.67) 

t  =  2.5  his 


6)   Next,  the  diving-gas  properties  were  determined. 
Since  the  diving-gas  was  air,  the  molecular  weights  and 
specific  heats  were  29.0  lb/ lb-mole  and  6.95  BTU/ lb-mole 
respectively.    Density  was  calculated  from  the  ideal 
equation  of  state,  equation  2.31. 
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=  [(29.0)  (36.9)  (144)] 
P      [  (1544)  (554.67) ] 

p  =  0.18002  lb/ ft3 


7)  The  psychrometric  properties  were  calculated  from 
equations  2.39  -  2.48.  The  psychrometric  chart  created 
for  this  case  is  presented  as  Figure  4.5. 

Ps  =  0.81564  psia  (from  Table  2-7) 
Since  the  operating  pressure  was  less  than  100  psia,  no 
correction   was   required   for   the   saturation   vapor 
pressure. 


x      =  0.81564 

ws  "    36  .9  -  0.81564 

Xws   =  0.0226  moles  water  vapor  /mole  dry  gas 


(X^  is  shown  as  point  A  on  Figure  4.5) 

Xwl   =    (0.95)  (0.0226) 

Xwl  =   0.0215  moles  water  vapor/ mole  dry  gas 


(X^  is  shown  as  point  B  on  Figure  4.5) 
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Moisture  Mole  Fraction 
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Xw2  =    (0.50)  (0.0226) 

Xw2  =  0.0113  moles  water  vapor/mole  dry  gas 

(X^  is  shown  as  point  D  on  Figure  4.5) 


H  =    (20,000) (0.0226) 

H  =  452.0  BTU/mole  dry  gas 


Hx   =    (0.95)  (452.0) 

Hx  =  429.4  BTU/mole  dry  gas 


H2   =    (0.50)  (452.0) 

H2  =   226  .0  BTU/mole  dry  gas 


8)  The  maximum  coil  temperature  was  determined  from  the 
psychrometric  chart  as  the  intersection  of  X^  and  the 
100%  saturation  curve.  This  point  is  at  approximately 
7  3.7°  F  and  is  shown  as  point  C  on  Figure  4.5. 

Tc  =  73.7°  F 

9)  The  energy  requirements  were  calculated  from 
equations  2.51  -  2.56. 
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H     =     [(6.95) (95)    +  429.4] 
a  29.0 

H.  =  37.6  BTU/lb 


H     =     [  (6.95)  (73.7)    +  226 .0] 
c  '  29.0 

Hr  =  25.5  BTU/lb 


Qd  =    (0.18002)  (39.9)  (37.6    -  25.5 
Qd  =  86.9  BTU/ min 


Ah  =     [6.95) (80   -  73.7)] 
29.0 

Ah  =  1.51  BTU/lb 


Qz  =    (0.18002)  (39.9)  (1.51 
Qr  =  10.8  BTU  I 'min 


0t  =  86.9   +10.8 
Qt  =  97  .7  £r£//min 


10)         The    moisture     removal     rate    was     calculated     from 
equations   2.57   -   2.59. 
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W,    =    (0.18002)  (  18.0158) 
1  29  .  0 

W1  =  0.013  36   lbjlbd 


W2   =    (0.18002)  (  18-0158) 
2  29.0 

W2  =  0.007  02  lbj lbd 


MRR  =    (0.01336    -  0.00702)  (39.9)  (0.18002)  (1440) 
MRR  =  65.5  lb  I day 


The  actual  program  printout  is  provided  as  Table  4-6.  The 
slight  discrepancies  between  the  example  and  the  actual  output 
were  caused  by  round-off  error  obtained  while  working  the 
example. 
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Table  4-6   Output  for  Underwater  Habitat 


UNDERWATER  HABITAT  SIMULATION 
Input  Conditions: 

Pressure  (psia) :  36.9 

Initial/Desired  Temperature  (°  F)  :  95.0/80.0 

Initial/Desired  Relative  Humidity  (%) :  95.0/50.0 

#  of  Divers/Working  Conditions:  4.0/3 

Chamber  Volume  (ft3):  2000.0 

Output  Conditions: 

Diving-Gas  Mixture:  AIR 

Free  Ventilation  Rate  at  Depth  (cfm) :  39.9 

Free  Ventilation  Rate  (SCFM) :  95.8 

Time  to  95%  Steady-State  (hrs) :  2.5 

Molecular  Weight  ( lb/ lb-mole) :  29.0 

Gas  Mixture  Density  (lb/ft3):  0.18002 

Maximum  Coil  Temperature  (°  F)  :  73.8 

Dehumidif ication  Energy  (BTU/min) :  87.0 

Reheat  Energy  (BTU/min):  10.8 

Total  Energy  Requirement  (BTU/min):  97.8 

Moisture  Removal  Rate  (lb/day):  65.4 
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CHAPTER  V 
SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A  FORTRAN  computer  program  has  been  developed  to  simulate 
a  hyperbaric  chamber  environment.  Although  hyperbaric 
chambers  have  been  used  for  many  years,  increasing  complexity 
requires  innovative  design  procedures.  This  program  allows 
the  user  to  design  life  support  systems,  temperature  and 
humidity  control  equipment,  and  determine  energy  requirements 
for  hyperbaric  chambers,  submarine  diver  lockout  chambers,  and 
underwater  habitats.  By  using  this  program,  the  design 
engineer  can  simulate  a  variety  of  scenarios  and  determine  the 
design  criteria  for  a  specific  chamber. 

Conclusions: 

1)  Acceptable  diving-gas  mixtures  are  obtained  by  ensuring 
the  oxygen  partial  pressure  of  the  mixture  remains  between 
0.16  and  1.6  atm  for  short  periods  of  time  and  below  0.6  atm 
for  extended  operations. 

2)  In  order  to  eliminate  the  narcosis  effects  of  nitrogen 
on  the  diver,  helium-oxygen  mixtures  are  used  below  150  fsw. 

3)  The  partial  pressure  of  carbon  dioxide  (C02)  must  be 
maintained  below  0.005  atm  during  extended  missions  in  order 
to  eliminate  discomfort  and  possible  injury  to  the  divers. 
This  can  be  accomplished  by  ventilation  of  the  chamber  or  by 
using  a  C02  absorption  system.  Ventilation  is  only 
recommended  when  air  is  used  as  the  diving-gas  due  to  the  cost 
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of  helium. 

4)  A  C02  absorption  system  consists  of  canisters  filled 
with  sodasorb.  The  diving-gas  mixture  flows  through  these 
canisters  and  the  C02  is  removed.  Although  this  program  can 
be  used  to  estimate  the  duration  of  the  C02  canisters  if  the 
amount  of  sodasorb  contained  is  known,  a  chemical  detection 
system  should  always  be  used  to  determine  when  to  replace  the 
absorbent  material. 

5)  The  molecular  weight  and  specific  heat  at  constant 
pressure  of  the  diving-gas  mixture  can  be  determined  from  the 
component  mole  fractions. 

6)  The  density  of  helium-oxygen  mixtures  is  determined  from 
the  virial  equation  of  state.  Three  virial  coefficients  are 
required  to  accurately  calculate  density.  Although,  the  first 
two  virial  coefficients  have  been  accurately  determined  for 
helium  and  oxygen,  only  one  second  virial  coefficient  for 
helium-oxygen  mixtures  has  been  determined.  Additionally,  no 
third  virial  coefficients  for  helium-oxygen  mixtures  have  been 
determined  to  date.  Therefore,  extremely  accurate  density 
calculations  for  helium-oxygen  mixtures  are  not  obtainable  at 
this  time. 

7)  Psychrometric  charts  have  been  used  for  years  to 
determine  dehumidif ication  and  cooling  requirements. 
Relationships  were  developed  to  compute  moisture  mole 
fractions  (X^)  allowing  a  computer  simulation  without  using 
the  psychrometric  charts.  This  program  is  easier  to  use  and 
more  accurate  than  psychrometric  charts. 
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8)  Enthalpy  of  moisture  content  (H)  is  easily  calculated 
once  X^  is  determined.  Enthalpy  is  an  important 
thermodynamic  property  for  calculating  energy  requirements. 

9)  Dehumidif ication  and  energy  requirements  can  be 
determined  for  varying  conditions  if  the  initial  and  desired 
temperatures  and  relative  humidities  are  known.  However,  the 
total  energy  requirement  neglects  heat  loss  due  to  heat 
transfer  between  the  chamber  and  the  environment. 
Additionally,  the  coil  temperature  required  to  perform  the 
desired  dehumidif ication  and  the  moisture  removal  rate  can  be 
determined. 

Recommendations : 

1)  Pressure  vessel  design  of  a  hyperbaric  chamber  could  be 
incorporated  into  a  computer  program.  NAVFAC  DM-39  (1982) 
provides  guidance  for  pressure  vessel  design  of  hyperbaric 
chambers  and  the  Pressure  Vessel  Handbook  (1977)  simplifies 
the  calculations  of  pressure  vessel  design  for  desired 
conditions. 

2)  Fire  protection  and  extinguishing  systems  along  with 
other  safety  design  criteria  could  be  incorporated  into  a 
computer  program.  Additionally,  electrical,  utility,  and 
communication  systems  could  be  considered.  NAVFAC  DM-39 
(1982)  provides  guidance  for  these  systems. 

3)  Different  diving-gas  mixtures  could  be  investigated  with 
additional  facility  requirements  incorporated  into  the 
program.     The   NOAA   Diving   Manual   (1991)   lists   some 
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alternatives  to  helium  and  research  into  the  effects  and 
advantages  of  these  gases  is  ongoing. 

4)  In  order  to  accurately  calculate  the  density  of  diving- 
gas  mixtures,  the  second  and  third  virial  coefficients  need  to 
be  accurately  determined.  Although  these  coefficients  exist 
for  helium  and  oxygen,  they  have  not  been  accurately 
determined  for  the  mixtures.  Texas  A  &  M  University  has  the 
facilities  to  accurately  determine  these  coefficients  if 
funding  is  available. 

5)  Heat  transfer  between  the  hyperbaric  chamber  and  the 
environment  is  not  considered  in  this  program.  A  detailed 
study  of  materials  and  insulation  would  be  necessary  prior  to 
incorporating  it  into  a  computer  program. 
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APPENDIX  A 
NOMENCLATURE 

B  -  second  virial  coefficient,  ft3/mole 

C  -  third  virial  coefficient,  ft6/mole2 

C02p  -  carbon  dioxide  production,  SCFM 

Cp  -  specific  heat  at  constant  pressure,  BTU/ lb-mole 

D  -  fourth  virial  coefficient,  ft9/mole3 

d  -  working  depth,  feet  of  seawater  (fsw) 

dmax  -  maximum  allowable  depth  for  a  diving-gas  mixture,  fsw 

dmin  -  minimum  allowable  depth  for  a  diving-gas  mixture,  fsw 

H  -  molar  enthalpy  of  moisture  content  at  saturation, 

BTU/mole  dry  gas 
Ha  -  moist-gas-mixture  enthalpy  at  initial  temperature,  BTU/ lb 
Hc  -  moist-gas-mixture  enthalpy  at  coil  temperature,  BTU/ lb 
H-l  -  molar  enthalpy  of  moisture  content  at  initial  conditions, 

BTU/mole  dry  gas 
H2  -  molar  enthalpy  of  moisture  content  at  desired  conditions, 

BTU/mole  dry  gas 
lbco2  -  carbon  dioxide  produced,  lb 
lbg  -  sodasorb  required  for  absorption,  lb 
m  -  mass  flowrate  of  carbon  dioxide,  lb/min 
MRR  -  moisture  removal  rate,  lb/day 

MW  -  molecular  weight  of  diving-gas  mixture,  lb/ lb-mole 
n  -  number  of  moles 
Nd  -  number  of  divers 
02c  -  oxygen  consumption,  SCFM 
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02r  -  oxygen  required  to  sustain  divers,  SCFM 

P  -  absolute  pressure,  psia 

Pmax  -  maximum  allowable  pressure  for  a  diving-gas  mixture, 

psia 
Pmin  -  minimum  allowable  pressure  for  a  diving-gas  mixture, 

psia 
PP02  -  partial  pressure  of  oxygen,  atm 
PC02  -  partial  pressure  of  carbon  dioxide,  atm 
Pg  -  vapor  pressure  of  water  in  saturated  moist  air,  psia 
%Heideal  -  percentage  of  helium  in  ideal  diving-gas  mixture 
%°2ideal  ~  percentage  of  oxygen  in  ideal  diving-gas  mixture 
%02max  -  maximum  acceptable  oxygen  percentage 
%02min  -  minimum  acceptable  oxygen  percentage 
Qd  -  energy  required  for  dehumidif ication,  BTU/min 
Qr  -  energy  required  to  reheat  from  coil  to  desired 

temperature,  BTU/min 
Qt  -  total  energy  required,  BTU/min 
R  -  gas  constant  for  C02 ,  35.11  lb/lbmole-R 
RHd  -  desired  relative  humidity,  % 
RH^   -  initial  relative  humidity,  % 
Rq  -  respiratory  quotient 

R^,  -  universal  gas  constant,  1544  ft-lb/lbmole-R 
Sair  -  ventilation  rate  of  air  at  standard  conditions,  SCFM 
T  -  initial  temperature,  °  F 

t  -  time  to  reach  95%  steady-state  conditions,  hr 
Tc  -  maximum  coil  temperature  for  dehumidif ication,  °  F 
Td  -  desired  temperature,  °  F 
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td  -  mission  duration,  hr 
U  -  internal  energy,  BTU 

V  -  volume,  ft3 

Vair  -  volumetric  flowrate  at  depth,  ft3/min 
Vch  -  volume  of  the  chamber,  ft3 

VI  -  specific  volume,  ft3/ lb 
Vm  -  molar  volume,  ft3/mole 

W2  -  humidity  ratio  at  initial  relative  humidity,  lbw/lbd 

W2  -  humidity  ratio  at  desired  relative  humidity,  lbw/lbw 

XHe  -  mole  fraction  of  helium 

X02  -  mole  fraction  of  oxygen 

X^  -  mole  fraction  of  water  at  saturation, 

moles  water  vapor/mole  dry  gas 
Xwl  -  moisture  mole  fraction  at  initial  conditions, 

moles  water  vapor/mole  dry  gas 
X^  ~  moisture  mole  fraction  at  desired  conditions, 

moles  water  vapor/mole  dry  gas 
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APPENDIX  B 
SOURCE  CODE 
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C  MAJOR  REPORT  -  "A  COMPUTER  SIMULATION  OF  A  HYPERBARIC 
C  CHAMBER  ENVIRONMENT." 

REAL   VAR,  DEPTH,  PRESS,  TFI ,  TRI ,  TFD,  RHI ,  RHD, 
+NUMDIV,  CHVOL,  FLOWRATE,  DUR,  VAR1 ,  02CONS,  PER02 , 
+PERHE,  02MIN,  02MAX,  DMIN,  DMAX,  PSIAMIN,  PSIAMAX,  VAIR, 
+SURAIR,  02REQ,  TIME,  MW,  CP,  DENS,  BC,  BF,  XWS,  XW1, 
+XW2,  H,  HI,  H2,  HA,  HC,  CT,  QD,  QR,  QT,  MRR,  AIR,  LBC02 , 
+LBS,  CASES 

C  VAR  =  A  VARIABLE  TO  DETERMINE  WHETHER  THE  USER  WISHES  TO 

C  CALCULATE  THE  GAS  MIXTURE  BASED  UPON  DEPTH  OR 

C  PRESSURE. 

C  DEPTH  =  A  WORKING  DEPTH  USED  TO  CALCULATE  THE  DIVING-GAS 

C  MIXTURE  IN  FEET. 

C  PRESS  =  PRESSURE  AT  DEPTH  IN  PSIA. 

C  TFI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 

C  FARENHEIT. 

C  TRI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 

C  RANKINE. 

C  TFD  =  THE  DESIRED  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 

C  FARENHEIT . 

C  RHI  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER. 

C  RHD  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER. 

C  NUMDIV  =  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

C  CHVOL  =  THE  INTERNAL  VOLUME  OF  THE  CHAMBER  IN  CUBIC 

C  FEET. 

C  FLOWRATE  =  THE  FLOWRATE  INSIDE  THE  CHAMBER  USED  FOR 

C  RECIRCULATING  THE  DIVING-GAS  MIXTURE  WHEN  C02 

C  SCRUBBERS  ARE  USED. 

C  DUR  =  THE  MISSION  DURATION  IN  HOURS. 

C  VAR1  =  A  VARIABLE  USED  TO  DETERMINE  THE  WORKING 

C  CONDITIONS. 

C  02CONS  =  OXYGEN  CONSUMPTION  IN  CUBIC  FEET  PER  MINUTE. 

C  PER02  =  THE  PERCENT  OF  OXYGEN  IN  THE  DIVING-GAS  MIXTURE. 

C  PERHE  =  THE  PERCENT  OF  HELIUM  IN  THE  DIVING-GAS  MIXTURE. 

C  02MIN  =  MINIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 

C  DIVING-GAS . 

C  02MAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 

C  DIVING-GAS. 

C  DMIN  =  THE  MINIMUM  DEPTH  IN  FEET  FOR  WHICH  A  GIVEN 

C  DIVING-GAS  CAN  BE  USED. 

C  DMAX  =  THE  MAXIMUM  DEPTH  IN  FEET  FOR  WHICH  A  GIVEN 

C  DIVING-GAS  CAN  BE  USED. 

C  PSIAMIN  =  THE  MINIMUM  PRESSURE  IN  PSIA  FOR  WHICH  A  GIVEN 

C  DIVING-GAS  CAN  BE  USED. 

C  PSIAMAX  =  THE  MAXIMUM  PRESSURE  IN  PSIA  FOR  WHICH  A  GIVEN 

C  DIVING-GAS  CAN  BE  USED. 

C  VAIR  =  VENTILATION  RATE  AT  DEPTH  IN  CUBIC  FEET  PER 

C  MINUTE. 

C  SURAIR  =  VENTILATION  RATE  AT  THE  SURFACE  IN  STANDARD 

C  CUBIC  FEET  PER  MINUTE. 
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C      02REQ  =  OXYGEN  REQUIRED  TO  SUSTAIN  DIVERS  IN  STANDARD 

C  CUBIC  FEET  PER  MINUTE  IF  C02  ABSORBERS  ARE  USED. 

C  TIME  =  THE  TIME  REQUIRED  FOR  THE  CHAMBER  TO  REACH  99% 

C  STEADY-STATE  CONDITIONS. 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 

C  (LB /LB-MOLE) . 

C      CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 

C  (BTU/MOLE) . 

C  DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU  FT) . 

C      BC  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  CM/G-MOLE) . 

C  BF  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  FT/LB-MOLE) . 

C  XWS  =  THE  MOISTURE  MOLE  FRACTION  OF  WATER  VAPOR  AT 

C  SATURATION  (MOLES  WATER  VAPOR/MOLE  DRY  GAS) . 

C  XW1  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 

C  HUMIDITY. 

C  XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 

C  HUMIDITY. 

C  H  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  SATURATION 

C  (BTU/MOLE  DRY  GAS) . 

C  HI  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 

C  INITIAL  R.H. 

C  H2  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 

C  DESIRED  R.H. 

C  HA  =  THE  MO I ST -GAS -MIXTURE  ENTHALPY  AT  INITIAL 

C  CONDITIONS  (BTU/LB  DRY  GAS) . 

C  HC  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  DESIRED 

C  CONDITIONS. 

C  CT  =  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION 

C  IN  DEGREES  FARENHEIT. 

C  QD  =  THE  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  (BTU/MIN) . 

C  QR  =  THE  ENERGY  REQUIRED  FOR  REHEAT  (BTU/MIN) . 

C  QT  =  THE  TOTAL  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  AND 

C  REHEAT  (BTU/MIN) . 

C  MRR  =  THE  REQUIRED  MOISTURE  REMOVAL  RATE  FOR  THE 

C  DEHUMIDIFICATION  PROCESS  (LB/ DAY) . 

C  LBC02  =  AMOUNT  OF  C02  PRODUCED  (LBS) . 

C  AIR  =  A  VARIABLE  TO  DETERMINE  WHETHER  OR  NOT  AIR  IS 

C  BEING  USED  AS  THE  DIVING-GAS. 

C  LBS  =  AMOUNT  OF  SODASORB  REQUIRED  FOR  C02  ABSORPTION 

C  (LBS)  . 

C  CASES  =  THE  NUMBER  OF  CASES  TO  BE  CALCULATED. 


C     THIS  PORTION  OF  THE  PROGRAM  WILL  ALLOW  THE  USER  TO  INPUT 
C      THE  NECESSARY  PARAMETERS  FROM  A  DATA  FILE. 

OPEN  (UNIT=3,FILE='VAR.DAT' , STATUS= ■ OLD ' ) 

OPEN  (UNIT=4 , FILE= ' OUTPUT . DAT ' , STATUS= ■ UNKNOWN ' ) 


READ  (3,*)  CASES 
DO  500  I  =  1, CASES 
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READ  (3,*)  VAR,  DEPTH,  TFI,  TFD,  RHI ,  RHD 
READ  (3,*)  NUMDIV,  CHVOL,  FLOWRATE,  DUR,  VAR1 


IF  (VAR.EQ.O)  THEN 

CALL  CONVERT (DEPTH,  PRESS) 

GO  TO  10 
ELSE 

PRESS=DEPTH 

GO  TO  10 
END  IF 
10     TRI=TFI+459.67 

IF  (VAR1.EQ.1)  THEN 

O2CONS=0.5/28.3 
ELSEIF  (VAR1.EQ.2)  THEN 

O2CONS=0.7/28.3 
ELSEIF  (VAR1.EQ.3)  THEN 

02CONS=1.5/28.3 
ELSEIF  (VAR1.EQ.4)  THEN 

O2CONS=2.0/28.3 
ELSE 

O2CONS=3.0/28.5 
END  IF 


C      THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  IDEAL  AND 
C      ACCEPTABLE  DIVING-GAS  MIXTURES. 

CALL  GASMIX  (PRESS, 02MIN, 02MAX) 

IF  (O2MAX.GE.21.0)  THEN 
AIR=1 
PERO2=21.0 

ELSE 

AIR=0 

IF  (O2MAX.GE.15.0)  THEN 

PERO2=15.0 
ELSEIF  (O2MAX.GE.10.0)  THEN 

PERO2=10.0 
ELSEIF  (O2MAX.GE.5.0)  THEN 

PERO2=5.0 
ELSEIF  (O2MAX.GE.3.0)  THEN 
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PER02=3 . 0 
ELSEIF  (02MAX.GE.1.0)  THEN 

PER02=1.0 
ELSEIF  (O2MAX.GE.0.5)  THEN 

PERO2=0.5 
ELSE 

PER02=0.1 
END  IF 

PERHE=100 . 0-PER02 

CALL  GASRANGE (PER02 , DMIN, DMAX, PSIAMIN, PSIAMAX) 

ENDIF 

C      THIS  PORTION  OF  THE  PROGRAM  WILL  CALCULATE  THE  REQUIRED 
C      VENTILATION  RATE  WITH  OR  WITHOUT  CARBON  DIOXIDE 
C      ABSORBERS . 

IF  (AIR.EQ.l)  THEN 

CALL  FREERATE (02 CONS , NUMDIV, TRI , PRESS , VAIR, SURAIR) 
ELSE 

CALL  ABSRATE (NUMDIV, 02REQ, 02 CONS , DUR, LBC02 , LBS) 

ENDIF 

CALL  STSTATE ( CHVOL , 02 CONS , NUMDIV , TRI , TIME ) 

C      THIS  PORTION  OF  THIS  PROGRAM  CALCULATES  THE  PROPERTIES 
C      OF  THE  DIVING-GAS  MIXTURE. 

CALL  M0LWT(PER02,MW,AIR) 

CALL  SPECHEAT(PER02,CP,AIR) 

CALL  DENSITY ( PER02 , TRI , PRESS , MW , DENS , BC , BF , AIR) 

C     THIS  PORTION  OF  THE  PROGRAM  DETERMINES  THE  PSYCHROMETRIC 
C      PROPERTIES. 

CALL  PSYCHART(TFI, TRI, PRESS, XWS,H) 

CALL  XWSH (RHI , RHD , XWS , XW1 , XW2 , H , HI , H2 ) 

CALL  C0ILTEMP(XW2, PRESS, CT) 

CALL  MGMENTH ( TFI , CT , CP , HI , H2 , MW , HA , HC ) 
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C      THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  ENERGY 

C     REQUIREMENTS  FOR  DEHUMIDIFICATION  AND  REHEAT  ALONG  WITH 

C      MOISTURE  REMOVAL  REQUIREMENTS. 

CALL  ENERGY ( DENS , AIR , VAIR , FLOWRATE , HA , HC , TFD , CT , CP , MW 
+QD,QR,QT) 

CALL  MOI STREM ( XW1 , XW2 , AIR , FLOWRATE , MW , VAIR , DENS ,  MRR) 

C      THIS  PORTION  OF  THE  PROGRAM  PRODUCES  THE  OUTPUT. 

IF  (AIR.EQ.O)  THEN 

WRITE  (4,12) 
12    FORMAT  (IX,'******************************************* 
+****************************  i  ) 

WRITE  (4,15)  I 
15     FORMAT  (35X, 'CASE  ',12) 

WRITE  (4,20) 
20    FORMAT  (IX,'******************************************* 
+****************************• ) 

WRITE  (4,25) 
25     FORMAT  ( 2 5X, ' "HYPERBARIC  CHAMBER  SIMULATION" » ) 

WRITE  (4,30) 
30    FORMAT  (IX,'******************************************* 

WRITE  (4,35) 
3  5     FORMAT  (5X,  ' INPUT  CONDITIONS:') 

WRITE  (4,40)  PRESS 
40     FORMAT  (5X, 'PRESSURE  (PSIA) : ' , 24X, 1F6 . 1) 

WRITE  (4,45)  TFI,TFD 
45    FORMAT  (5X,  ■  INITIAL/DESIRED  TEMP  (DEG  F)  :  '  ,  9X,  1F4  . 1,  '  /  '  , 
+1F4.1) 

WRITE  (4,50)  RHI,RHD 
50    FORMAT  (5X,  ' INITIAL/DESIRED  REL  HUM  (%):', 10X, 1F4 . 1 ,'/' , 
+1F4.1) 

WRITE  (4,55)  NUMDIV,VAR1 
55     FORMAT  (5X,'#  OF  DIVERS/WORKING  CONDITIONS :', 7X, 1F4 . 1 , 
+'/ ' ,1F2.0) 

WRITE  (4,60)  CHVOL 
60     FORMAT  ( 5 X, ■ CHAMBER  VOLUME  (CU  FT) : ' , 17X, 1F6 . 1) 
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WRITE  (4,70)  FLOWRATE 
70     FORMAT  (5X,'FLOW  RATE  (CU  FT/MIN) : • , 18X, 1F6 . 1) 

WRITE  (4,75)  DUR 
75     FORMAT  (5X, 'MISSION  DURATION  (HRS) : • , 17X, 1F6 . 1) 

WRITE  (4,100) 
100   FORMAT  (IX,'******************************************* 

WRITE  (4,105) 
105    FORMAT  (5X, 'OUTPUT: ' ) 

WRITE  (4,110)  PERHE,PER02 
110    FORMAT  (5X, 'DIVING-GAS  MIXTURE  (%HE/%02 ) : ' , 9X, 1F4 . 1 , 
+'/ ' ,1F4.1) 

WRITE  (4,115)  DMIN,DMAX 
115    FORMAT  (5X, 'DIVING-GAS  MIXTURE  RANGE  (FSW) : ' , 7X, 1F6 . 1, 
+  '/ ' ,1F6.1) 

WRITE  (4,120)  02REQ 
120    FORMAT  (5X, 'OXYGEN  REQUIRED  (SCFM) : ' , 18X, 1F8 . 4 ) 

WRITE  (4,125)  LBC02 
125    FORMAT  (5X,'C02  PRODUCED  (LB) : ' , 22X, 1F6 . 1) 

WRITE  (4,130)  LBS 
130    FORMAT  (5X, ' SODASORB  REQUIRED: ', 22X, 1F6 . 1) 

WRITE  (4,135)  TIME 
135    FORMAT  (5X, 'TIME  TO  95%  STEADY  STATE  (HRS) : ' , 9X, 1F6 . 1) 

WRITE  (4,140)  MW 
140    FORMAT  ( 5 X, 'MOLECULAR  WEIGHT  (LB/MOLE) :', 13X, 1F8 . 3 ) 

WRITE  (4,145)  DENS 
145    FORMAT  (5X,'GAS  MIXTURE  DENSITY  (LB/CU  FT) : ' , 9X, 1F10 . 5) 

WRITE  (4,150)  CT 
150   FORMAT  (5X,  'MAXIMUM  COIL  TEMPERATURE  (DEG  F)  :  ',7X,1F6.1) 

WRITE  (4,155)  QD 
155    FORMAT  (5X, ' DEHUMIDIFICATION  ENERGY  (BTU/MIN) : ' , 6X, 
+1F6.1) 

WRITE  (4,160)  QR 
160    FORMAT  (5X, 'REHEAT  ENERGY  (BTU/MIN) :', 16X, 1F6 . 1) 

WRITE  (4,165)  QT 
165    FORMAT  (5X, ' ENERGY  REQUIREMENTS  (BTU/MIN) :', 10X, 1F6 . 1) 
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WRITE  (4,170)  MRR 
170    FORMAT  (5X, 'MOISTURE  REMOVAL  RATE  (LB/DAY) :', 9X, 1F6 . 1) 

WRITE  (4,175) 
175   FORMAT  (IX,1******************************************* 

ELSE 

WRITE  (4,200) 
2  00   FORMAT  (IX,'******************************************* 
+**************************** • ) 

WRITE  (4,205) 
205    FORMAT  ( 2 5X, ' "HYPERBARIC  CHAMBER  SIMULATION" ' ) 

- 

WRITE  (4,210) 
210   FORMAT  (IX,'******************************************* 

WRITE  (4,215) 
215    FORMAT  (5X, ' INPUT  CONDITIONS: ' ) 

WRITE  (4,220)  PRESS 
220    FORMAT  (5X, • PRESSURE  (PSIA) : ' , 24X, 1F6 . 1) 

WRITE  (4,225)  TFI,TFD 
225   FORMAT  (5X,  '  INITIAL/DESIRED  TEMP  (DEG  F)  :  •  ,  9X,  1F4  .  1 ,  •  /  '  , 
+1F4.1) 

WRITE  (4,230)  RHI,RHD 
230   FORMAT  (5X,  • INITIAL/DESIRED  REL  HUM  (%):', 10X, 1F4 . 1,  »/' , 
+1F4.1) 

WRITE  (4,235)  NUMDIV,VAR1 
235    FORMAT  (5X,'#  OF  DIVERS/WORKING  CONDITIONS :', 7X, 1F4 . 1 , 
+•/' ,1F2.0) 

WRITE  (4,240)  CHVOL 
240    FORMAT  ( 5 X, • CHAMBER  VOLUME  (CU  FT) : • , 17X, 1F6 . 1) 

WRITE  (4,245) 
245   FORMAT  (IX,'******************************************* 

WRITE  (4,250) 
250    FORMAT  (5X, ' OUTPUT: • ) 

WRITE  (4,255) 
255    FORMAT  (5X, 'DIVING-GAS  MIXTURE: ', 23X, 'AIR' ) 
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WRITE  (4,260) 
260    FORMAT  (5X, ' DIVING-GAS  MIXTURE  RANGE  (FSW):',10X, 
+'0.0/60.0' ) 

WRITE  (4,265)  VAIR 

2  65    FORMAT  (5X,  'FREE  VENTILATION  RATE  AT  DEPTH  (CFM)  :  '  ,3X, 

+1F6.1) 

WRITE  (4,270)  SURAIR 
270    FORMAT  (5X, ' FREE  VENTILATION  RATE  (SCFM) : ' , 11X, 1F6 . 1) 

WRITE  (4,275)  TIME 
275    FORMAT  (5X,'TIME  TO  95%  STEADY-STATE  (HRS) : » , 9X, 1F6 . 1) 

WRITE  (4,280)  MW 
280    FORMAT  (5X, 'MOLECULAR  WEIGHT  (LB/MOLE) :', 15X, 1F4 . 1) 

WRITE  (4,285)  DENS 
285    FORMAT  (5X, 'DENSITY  OF  AIR  (LB/CU  FT) : ' , 14X, 1F10 . 5) 

WRITE  (4,290)  CT 
290   FORMAT  (5X,  'MAXIMUM  COIL  TEMPERATURE  (DEG  F)  :  '  ,7X,1F6.1) 

WRITE  (4,295)  QD 
295    FORMAT  (5X, ' DEHUMIDIFICATION  ENERGY  (BTU/MIN) : ' , 6X, 
+1F6.1) 

WRITE  (4,300)  QR 

3  00    FORMAT  (5X, 'REHEAT  ENERGY  (BTU/MIN) :', 16X, 1F6 . 1) 

WRITE  (4,305)  QT 
305    FORMAT  (5X, ' ENERGY  REQUIREMENTS  (BTU/MIN) :', 10X, 1F6 . 1) 

WRITE  (4,310)  MRR 
310    FORMAT  (5X, 'MOISTURE  REMOVAL  RATE  (LB/DAY) :', 9X, 1F6 . 1) 

WRITE  (4,315) 
315   FORMAT  (IX,'******************************************* 

END  IF 
500    CONTINUE 
END 
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SUBROUTINE  CONVERT ( DEPTH , PRESS ) 

C      THIS  SUBROUTINE  CONVERTS  A  WORKING  DEPTH  IN  FEET  TO  A 
C      PRESSURE  IN  PSIA. 

REAL  DEPTH,  PRESS 

C      DEPTH  =  A  WORKING  DEPTH  IN  FEET. 

C      PRESS  =  A  CONVERTED  PRESSURE  IN  PSIA. 

PRESS= ( (DEPTH+33 . 06633) / 33 . 06633) *14 . 696 

END 


SUBROUTINE  GASMIX ( PRESS , 02MIN , 02MAX) 

C  THIS  SUBROUTINE  CALCULATES  AN  ACCEPTABLE  DIVING  GAS 

C  MIXTURE  BASED  UPON  A  MINIMUM  OXYGEN  PARTIAL  PRESSURE  OF 

C  0.16  ATM  AND  A  MAXIMUM  OXYGEN  PARTIAL  PRESSURE  OF  0 . 6 

C  ATM. 

REAL  PRESS,  PATM,  01MIN,  02MIN,  01MAX,  02MAX 

C  PRESS  =  THE  CHAMBER  PRESSURE  IN  PSIA. 

C  PATM  =  THE  CHAMBER  PRESSURE  IN  ATM. 

C  01MIN  =  MINIMUM  ACCEPTABLE  OXYGEN  FRACTION  IN  THE 

C  DIVING-GAS. 

C  02MIN  =  MINIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 

C  DIVING-GAS. 

C  01MAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  FRACTION  IN  THE 

C  DIVING-GAS. 

C  02MAX  =  MAXIMUM  ACCEPTABLE  OXYGEN  PERCENTAGE  IN  THE 

C  DIVING-GAS . 

PATM=PRESS/ 14 . 696 
01MIN=0.16/PATM 
02MIN=01MIN*100 . 0 
01MAX= 0.6 /PATM 
O2MAX=OlMAX*100 . 0 

END 
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SUBROUTINE  GASRANGE(PER02 , DMIN, DMAX, PSIAMIN, PSIAMAX) 

C      THIS  SUBROUTINE  CALCULATES  A  DEPTH  AND  PRESSURE  RANGE 
C      FOR  A  GIVEN  DIVING-GAS  MIXTURE. 

REAL  PER02,  PMIN,  DMIN,  PMAX,  DMAX,  PSIAMIN,  PSIAMAX 

C  PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 

C  MIXTURE. 

C  PMIN  =  THE  MINIMUM  PRESSURE  (ATM)  BASED  UPON  AN  OXYGEN 

C  PARTIAL  PRESSURE  OF  0.16  ATM. 

C  DMIN  =  THE  MINIMUM  DEPTH  (FT)  FOR  A  GIVEN  DIVING-GAS. 

C  PMAX  =  THE  MAXIMUM  PRESSURE  (ATM)  BASED  UPON  AN  OXYGEN 

C  PARTIAL  PRESSURE  OF  0.6  ATM. 

C  DMAX  =  THE  MAXIMUM  DEPTH  (FT)  FOR  A  GIVEN  DIVING-GAS. 

C  PSIAMIN  =  THE  MINIMUM  PRESSURE  (PSIA)  FOR  A  GIVEN 

C  DIVING-GAS. 

C  PSIAMAX  =  THE  MAXIMUM  PRESSURE  (PSIA)  FOR  A  GIVEN 

C  DIVING-GAS. 

PMIN=16.0/PERO2 
DMIN=PMIN*33 . 06633-33 . 06633 
PMAX=60.0/PERO2 
DMAX=PMAX*33 . 06633-33 . 06633 
PSIAMIN=14 . 696*PMIN 
PSIAMAX=14 . 696*PMAX 

END 


SUBROUTINE  FREERATE  ( 02 CONS ,  NUMDIV ,  TRI ,  PRESS ,  VAIR ,  SURAIR) 

C     THIS  SUBROUTINE  CALCULATES  THE  REQUIRED  VENTILATION  RATE 
C      FOR  THE  CHAMBER  WITHOUT  CARBON  DIOXIDE  ABSORBERS. 

REAL  02 CONS,  NUMDIV,  TRI,  PRESS,  VAIR,  SURAIR,  C02PROD, 
+MFLOW 

C  02CONS  =  OXYGEN  CONSUMPTION  IN  CUBIC  FEET  PER  MINUTE. 

C  NUMDIV  =  THE  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

C  TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 

C  RANKINE. 

C  PRESS  =  THE  CHAMBER  PRESSURE  IN  PSIA. 

C  VAIR  =  THE  VENTILATION  RATE  AT  DEPTH  (CU  FT/MIN) . 

C  SURAIR  =  THE  VENTILATION  RATE  AT  THE  SURFACE  (SCFM) . 

C  C02PROD  =  THE  TOTAL  AMOUNT  OF  CARBON  DIOXIDE  PRODUCED 

C  (CU  FT/MIN) . 

C  MFLOW  =  THE  MASS  FLOW  RATE  OF  CARBON  DIOXIDE  (LB/MIN) . 
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C02PROD=0 . 9*02CONS*NUMDIV 
MFLOW=0. 11372474 *C02PR0D 
VAIR=3 . 31817426*MFLOW*TRI 
SURAIR=(VAIR*PRESS*529.67) / (14 . 696*TRI) 


END 


SUBROUTINE  ABSRATE ( NUMDI V , 02REQ , 02 CONS , DUR , LBC02 ,  LBS ) 

C     THIS  SUBROUTINE  CALCULATES  THE  REQUIRED  VENTILATION  RATE 
C     OF  THE  CHAMBER  IF  CARBON  DIOXIDE  ABSORBERS  ARE  UTILIZED. 

REAL  NUMDIV,  02REQ,  02 CONS,  DUR,  LBC02 ,  MFLOW,  LBS, 
+C02PR0D 

C  NUMDIV  =  THE  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

C  02REQ  =  THE  AMOUNT  OF  PURE  OXYGEN  REQUIRED  (SCFM) . 

C  02 CONS  =  THE  OXYGEN  CONSUMPTION  (SCFM) . 

C  DUR  =  THE  MISSION  DURATION  IN  HOURS. 

C  LBC02  =  THE  AMOUNT  OF  C02  PRODUCED  (LB) . 

C  MFLOW  =  THE  MASS  FLOWRATE  OF  CARBON  DIOXIDE  (LB/MIN) . 

C  C02PR0D  =  CARBON  DIOXIDE  PRODUCTION  (CU  FT/MIN) . 

C  LBS  =  AMOUNT  OF  SODASORB  REQUIRED  FOR  C02  ABSORPTION 

C  (LBS) . 

02REQ=NUMDIV*02CONS 
C02PROD=0 . 9*02CONS*NUMDIV 
MFLOW=0.11372474*CO2PROD 
LBCO2=MFLOW*60 . 0*DUR 
LBS=LBC02/ 0.245 


END 


SUBROUTINE  STSTATE ( CHVOL , 02 CONS , NUMDIV , TRI , TIME ) 

C      THIS  SUBROUTINE  CALCULATES  THE  TIME  REQUIRED  FOR  THE 
C      CHAMBER  TO  REACH  A  95%  STEADY-STATE  CONDITION. 

REAL  CHVOL,  02 CONS,  NUMDIV,  TRI,  TIME,  C02PR0D,  MFLOW, 
+VAIR 

C      CHVOL  =  THE  INTERIOR  CHAMBER  VOLUME  (CUBIC  FEET) . 

C      02C0NS  =  OXYGEN  CONSUMPTION  (CU  FT/MIN) . 

C      NUMDIV  =  THE  NUMBER  OF  DIVERS  IN  THE  CHAMBER. 

C  TRI  =  THE  INITIAL  TEMPERATURE  OF  THE  CHAMBER  IN  DEGREES 
C  RANKINE. 

C  TIME  =  THE  TIME  REQUIRED  TO  REACH  99%  STEADY-STATE 
C  (HRS) . 
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C  C02PR0D  =  THE  CARBON  DIOXIDE  PRODUCTION  (CU  FT/MIN) . 
C  MFLOW  =  THE  MASS  FLOWRATE  OF  CARBON  DIOXIDE  (LB/MIN) 
C      VAIR  =  THE  VENTILATION  RATE  AT  DEPTH  (CU  FT/MIN) . 

CO2PROD=0 . 9*02CONS*NUMDIV 
MFLOW=0.11372474*CO2PROD 
VAIR=3 . 31817426*MFLOW*TRI 

TIME=(-0. 0167 *CHVOL/VAIR)*LOG(1.0-(0. 2863 *VAIR) 
+   /(MFLOW*TRI) ) 

END 


SUBROUTINE  MOLWT (PER02 ,MW, AIR) 

C      THIS  SUBROUTINE  CALCULATES  THE  MOLECULAR  WEIGHT  OF  THE 
C      DIVING-GAS  MIXTURE. 

REAL  PER02,  MW,  X02  ,  XHE,  AIR 

C  PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 

C  MIXTURE. 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  (LB/ LB-MOLE) . 

C  X02  =  THE  MOLE  FRACTION  OF  OXYGEN  IN  THE  MIXTURE. 

C  XHE  =  THE  MOLE  FRACTION  OF  HELIUM  IN  THE  MIXTURE. 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

IF  (AIR.EQ.l)  THEN 

MW=29.0 

ELSE 

X02=PER02/ 100.0 
XHE=1.0-XO2 
MW=(XO2*32.00)+(XHE*4.002  6) 

ENDIF 

END 


SUBROUTINE  SPECHEAT ( PER02 , CP , AIR) 

C      THIS  SUBROUTINE  CALCULATES  THE  SPECIFIC  HEAT  OF  THE 
C      DIVING-GAS  MIXTURE. 

REAL  PER02,  CP,  X02 ,  XHE,  AIR 

C      PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 
C  MIXTURE. 
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C  CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 

C  (BTU/MOLE) . 

C  X02  =  THE  MOLE  FRACTION  OF  OXYGEN  IN  THE  MIXTURE. 

C  XHE  =  THE  MOLE  FRACTION  OF  HELIUM  IN  THE  MIXTURE. 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR, 

IF  (AIR.EQ.l)  THEN 

CP=6.95 

ELSE 

X02=PER02/ 100.0 

XHE=1.0-XO2 

CP=(X02*6.95)+(XHE*4.965) 

ENDIF 

END 


SUBROUTINE  DENSITY (PER02 , TRI , PRESS , MW , DENS , BC , BF , AIR) 

C     THIS  SUBROUTINE  CALCULATES  THE  DENSITY  OF  THE  DIVING-GAS 
C      MIXTURE. 

REAL  PER02,  TRI,  PRESS,  MW,  DENS,  BC,  BF,  CONST,  MOLVOL, 
+AIR 

C  PER02  =  THE  PERCENTAGE  OF  OXYGEN  IN  THE  DIVING-GAS 

C  MIXTURE. 

C  TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 

C  RANKINE. 

C  PRESS=  THE  INITIAL  CHAMBER  PRESSURE  IN  PSIA. 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 

C  (LB/ LB-MOLE) . 

C  DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU 

C  FEET) . 

C  BC  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  CM/MOLE) . 

C  BF  =  THE  SECOND  VIRIAL  COEFFICIENT  (CU  FT/MOLE) . 

C  CONST  =  RT/P 

C  MOLVOL  =  MOLAR  VOLUME  (CU  FT /MOLE) 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

IF  (AIR.EQ.l)  THEN 

DENS=(MW*PRESS*144.0) / (1544.0*TRI) 

ELSE 

X02=PER02/ 100.0 
XHE=1.0-XO2 
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BC=(11.44*XHE**2)-(8.8*XHE*X02)-(16.24*X02**2) 
BF=BC*3 . 5315E-05*454 . 0 
CONST=(1544.0*TRI) / (PRESS*144  .  0) 

MOLVOL=(CONST+SQRT( (CONST**2 ) + (4 . 0*CONST*BF) ) ) /2.0 
DENS=MW/MOLVOL 


ENDIF 
END 


SUBROUTINE  PS YCHART ( TFI , TRI , PRESS , XWS , H ) 

C      THIS  SUBROUTINE  CALCULATES  THE  MOISTURE  MOLE  FRACTION 
C      AND  ENTHALPY  OF  MOISTURE  UTILIZING  THE  PSYCHROMETRIC 
C      CHART. 

REAL  TFI,  TRI,  PRESS,  XWS,  H,  PS,  DPS,  VL 

C  TFI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 

C  FARENHEIT. 

C  TRI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 

C  RANKINE. 

C  PRESS  =  THE  INITIAL  CHAMBER  PRESSURE  IN  PSIA. 

C  XWS  =  THE  MOISTURE  MOLE  FRACTION  (MOLES  WATER  VAPOR/MOLE 

C  DRY  GAS) . 

C  H  =  THE  ENTHALPY  OF  MOISTURE  CONTENT  (BTU/MOLE  DRY  GAS) . 

C  PS  =  THE  SATURATION  PRESSURE  OF  WATER  (PSIA) . 

C  DPS  =  THE  CHANGE  IN  SATURATION  PRESSURE  OF  WATER  FOR 

C  PRESSURES  ABOVE  100  PSIA. 

C  VL  =  THE  SPECIFIC  VOLUME  OF  THE  SATURATED  LIQUID  (CU 

C  FT/LB) . 

C     THIS  PORTION  OF  THE  SUBROUTINE  CALCULATES  THE  SATURATION 
C      PRESSURE  OF  WATER  GIVEN  THE  DRY  BULB  TEMPERATURE. 

PS=(1.9285E-11*TFI**5)+(9.8683E-10*TFI**4)+ 
+(3.3193E-07*TFI**3)+(1.9916E-05*TFI**2)+ 
+(1.0488E-03*TFI)+2.2  080E-02 


C  THIS  PORTION  OF  THE  SUBROUTINE  CORRECTS  THE  SATURATION 
C  PRESSURE  FOR  PRESSURES  ABOVE  100  PSIA  USING  POYNTING • S 
C      EQUATION . 

IF  (PRESS. GT. 100.0)  THEN 

IF  (TFI. GT. 32.0. AND. TFI. LE. 52.0)  THEN 

VL=0. 01602 
ELSEIF  (TFI. GT. 52.0. AND. TFI. LE. 60.0)  THEN 

VL=0. 01603 
ELSEIF  (TFI.GT.60.0.AND.TFI.LE.66.0)  THEN 

VL=0. 01604 
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ELSEIF 

[TFI.GT. 

66. 

0 

.AND.TFI 

VL=0. 

01605 

ELSEIF 

(TFI.GT. 

71. 

0 

.AND.TFI 

VL=0. 

01606 

ELSEIF 

[TFI.GT. 

76. 

0 

.AND.TFI 

VL=0. 

01607 

ELSEIF 

[TFI.GT, 

80. 

0 

.AND.TFI 

VL=0. 

01608 

ELSEIF 

[TFI.GT. 

84. 

0 

.AND.TFI 

VL=0. 

01609 

ELSEIF 

[TFI.GT. 

87. 

0 

.AND.TFI 

VL=0. 

01610 

ELSEIF 

[TFI.GT. 

91. 

0 

.AND.TFI 

VL=0. 

01611 

ELSEIF  1 

[TFI.GT. 

94. 

0 

.AND.TFI 

VL=0, 

01612 

ELSEIF  1 

[TFI.GT. 

97. 

0 

.AND.TFI 

VL=0. 

01613 

ELSEIF  I 

[TFI.GT. 

100.0)  THEN 

VL=0. 

01614 

END  IF 

LE.71.0)  THEN 

LE.7  6.0)  THEN 

LE.80.0)  THEN 

LE.84.0)  THEN 

LE.87.0)  THEN 

LE.91.0)  THEN 

LE.94.0)  THEN 

LE.97.0)  THEN 
LE. 100.0)  THEN 


DPS=(1.67876*PS*VL*PRESS)/TRI 
PS=PS+DPS 

GO  TO  40 

ELSE 

GO  TO  40 

END  IF 

40     XWS=PS/ (PRESS-PS) 
H=20000.0*XWS 

END 


SUBROUTINE  XWSH (RHI , RHD , XWS , XW1 , XW2 , H , HI , H2 ) 

C      THIS  SUBROUTINE  CALCULATES  THE  INITIAL  AND  DESIRED 
C      MOISTURE  MOLE  FRACTIONS  AND  MOLAR  ENTHALPIES. 

REAL  RHI,  RHD,  XWS,  XW1,  XW2 ,  H,  HI,  H2 ,  RHI,  RH2 

C  RHI  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER  IN 

C  PERCENT . 

C  RHD  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER  IN 

C  PERCENT. 

C  XWS  =  THE  MOISTURE  MOLE  FRACTION  FOR  WATER  VAPOR  AT 

C  SATURATION. 


108 


C  XW1  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 

C  HUMIDITY. 

C  XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 

C  HUMIDITY. 

C  H  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT 

C  SATURATION. 

C  HI  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 

C  INITIAL  R.H. 

C  H2  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AN  THE 

C  DESIRED  R.H. 

C  RH1  =  THE  INITIAL  RELATIVE  HUMIDITY  IN  THE  CHAMBER 

C  (FRACTION) . 

C  RH2  =  THE  DESIRED  RELATIVE  HUMIDITY  IN  THE  CHAMBER 

C  (FRACTION) . 

RH1=RHI/ 100.0 

RH2=RHD/ 100.0 

XW1=RH1*XWS 

XW2=RH2*XWS 

H1=RH1*H 

H2=RH2*H 

END 


SUBROUTINE  COILTEMP (XW2 , PRESS, CT) 

C     THIS  SUBROUTINE  CALCULATES  THE  COIL  TEMPERATURE  REQUIRED 
C      FOR  DEHUMIDIFICATION  TO  A  DESIRED  RELATIVE  HUMIDITY. 


REAL  XW2,  PRESS,  PS,  CT 


C  XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 

C  HUMIDITY. 

C  PRESS  =  THE  INITIAL  CHAMBER  PRESSURE  (PSIA) . 

C  PS  =  THE  SATURATION  PRESSURE  OF  WATER  (PSIA) . 

C  CT  =  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION. 

PS= ( XW2  *PRESS ) / ( 1 . 0+XW2 ) 

CT=(2.4135E02*PS**5)-(8.4538E02*PS**4)+ 
+(1.1635E03*PS**3)-(8.1910E02*PS**2)+ 
+ (3 . 5407E02*PS) +7 . 0106 

END 
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SUBROUTINE  MGMENTH (TFI , CT , CP , HI , H2 , MW , HA , HC) 

C      THIS  SUBROUTINE  CALCULATES  THE  MOIST-GAS-MIXTURE 
C      ENTHALPIES  AT  THE  INITIAL  CONDITIONS  AND  COIL 
C      TEMPERATURE . 

REAL  TFI,  CT,  CP,  HI,  H2 ,  MW,  HA,  HC 

C  TFI  =  THE  INITIAL  CHAMBER  TEMPERATURE  IN  DEGREES 

C  FARENHEIT. 

C  CT  =  THE  REQUIRED  COIL  TEMPERATURE  FOR  DEHUMIDIFICATION 

C  (F). 

C  CP  =  THE  SPECIFIC  HEAT  OF  THE  GAS  MIXTURE  (BTU/MOLE) . 

C  HI  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 

C  INITIAL  R.H. 

C  H2  =  THE  MOLAR  ENTHALPY  OF  MOISTURE  CONTENT  AT  THE 

C  DESIRED  R.H. 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  (LB/LB-MOLE) . 

C  HA  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  INITIAL 

C  CONDITIONS. 

C  HC  =  THE  MOIST-GAS -MIXTURE  ENTHALPY  AT  THE  DESIRED 

C  CONDITIONS. 

HA= ( (CP*TFI) +H1) /MW 
HC= ( (CP*CT) +H2 ) /MW 

END 


SUBROUTINE  ENERGY ( DENS , AIR , VAIR , FLOWRATE , HA , HC , TFD , CT , 
+CP,MW,QD,QR,QT) 

C      THIS  SUBROUTINE  CALCULATES  THE  ENERGY  REQUIRED  FOR 
C      DEHUMIDIFICATION  AND  REHEAT. 

REAL  DENS,  AIR,  VAIR,  FLOWRATE,  HA,  HC,  TFD,  CT,  CP,  MW, 
+QD,  DH,  QR,  QT 

C  DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU  FT) . 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

C  VAIR  =  THE  VENTILATION  OR  CIRCULATION  RATE  AT  DEPTH  (CU 

C  FT/MIN) . 

C  FLOWRATE  =  THE  CIRCULATION  RATE  (CU  FT/MIN) . 

C  HA  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  INITIAL 

C  CONDITIONS. 

C  HC  =  THE  MOIST-GAS-MIXTURE  ENTHALPY  AT  THE  DESIRED 

C  CONDITIONS. 

C  TFD  =  THE  DESIRED  CHAMBER  TEMPERATURE  AFTER  REHEAT  (F) . 

C  CT  =  THE  COIL  TEMPERATURE  REQUIRED  FOR  DEHUMIDIFICATION 

110 


C  CP  =  THE  SPECIFIC  HEAT  OF  THE  DIVING-GAS  MIXTURE 

C  (BTU/MOLE) . 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  MIXTURE  (LB/LB-MOLE) . 

C  QD  =  THE  ENERGY  REQUIRED  FOR  DEHUMIDIFICATION  (BTU/MIN) . 

C  DH  =  THE  CHANGE  IN  ENTHALPY  FROM  THE  COIL  TEMPERATURE  TO 

C  THE  DESIRED  CHAMBER  TEMPERATURE  (BTU/MIN) . 

C  QR  =  THE  ENERGY  REQUIRED  FOR  REHEAT  (BTU/MIN) . 

C  QT  =  THE  TOTAL  ENERGY  REQUIRED  DEHUMIDIFICATION  AND 

C  REHEAT  (BTU/MIN) . 

IF  (AIR.EQ.O)  THEN 
VAIR=FLOWRATE 

ELSE 

VAIR=VAIR 

END  IF 

QD=DENS*VAIR* (HA-HC) 
DH=CP* (TFD-CT) /MW 
QR=DENS*VAIR*DH 
QT=QD+QR 


END 


SUBROUTINE  MOISTREM ( XW1 , XW2 , AIR , FLOWRATE , MW , VAIR , 
+DENS,MRR) 

C      THIS  SUBROUTINE  CALCULATES  THE  MOISTURE  REMOVAL  RATE 
C      REQUIRED  FOR  DEHUMIDIFICATION. 

REAL  Wl,  W2,  XW1,  XW2 ,  AIR,  FLOWRATE,  MW,  VAIR,  DENS, 
+MRR 

C  Wl  =  THE  HUMIDITY  RATIO  AT  THE  INITIAL  RELATIVE 

C  HUMIDITY. 

C  W2  =  THE  HUMIDITY  RATIO  AT  THE  DESIRED  RELATIVE 

C  HUMIDITY. 

C  XW1  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  INITIAL  RELATIVE 

C  HUMIDITY. 

C  XW2  =  THE  MOISTURE  MOLE  FRACTION  AT  THE  DESIRED  RELATIVE 

C  -  HUMIDITY. 

C  AIR  =  A  VARIABLE  TO  DETERMINE  IF  THE  DIVING-GAS  IS  AIR. 

C  FLOWRATE  =  THE  CIRCULATION  RATE  (CU  FT/MIN) . 

C  MW  =  THE  MOLECULAR  WEIGHT  OF  THE  DIVING-GAS  MIXTURE 

C  (LB/ LB-MOLE) . 

C  VAIR  =  THE  VENTILATION  RATE  AT  DEPTH  (CU  FT/MIN) . 

C  DENS  =  THE  DENSITY  OF  THE  DIVING-GAS  MIXTURE  (LB/CU  FT) . 

C  MRR  =  THE  REQUIRED  MOISTURE  REMOVAL  RATE  FOR  THE 

C  DEHUMIDIFICATION  PROCESS  (LB/DAY) . 


Ill 


IF  (AIR.EQ.O)  THEN 
VAIR=FLOWRATE 

ELSE 

VAIR=VAIR 

ENDIF 

W1=XW1* (18 . 0158 /MW) 
W2=XW2*(18.0158/MW) 
MRR=(W1-W2)*VAIR*DENS*1440.0 

END 
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